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ABSTRACT 



This study was undertaken to determine if the US Army's Ml 19 
Howitzer could be safely towed in the “firing" position using the 
existing A-Frame brace to support the gun tube during movement. A 
computer model of the howitzer was developed that would predict 
motions of the howitzer and consequently, loads on critical 

components. The results obtained from this computer model were 

validated by comparing them to actual data obtained from an 
instrumented test done on the howitzer at Aberdeen Proving 

Grounds, Md. The validated model was used to predict the "worst 
case" loads on the howitzer's A-Frame and the T-Bar during 

movement in the "firing" position. The maximum predicted stresses 
in the A-Frame and the T-Bar from the "worst case" loads were 
compared to the yield stress. These comparisons showed that the A- 
Frame and T-Bar were strong enough to withstand the "worst case" 
predicted loads that would be caused by towing the Ml 19 Howitzer 
in the "firing" position. 



m 



TABLE OF CONTENTS 



I. INTRODUCTION 1 

A. GENERAL 1 

B. BACKGROUND 1 

C. METHOD OF ANALYSIS 2 

II. HUMMV/HOWITZER SYSTEM COMPUTER MODEL 4 

A. HUMMV/HOWITZER SYSTEM PARAMETERS 4 

1. HUMMV and Howitzer Tires 7 

2. HUMMV Suspension System 7 

3. Howitzer Suspension System 8 

4. A-Frame/T-Bar Stiffness 17 

5 . Equilibrator Springs 17 

6. Elevating Clutch 17 

7. HUMMV Component Mass Properties 18 

8. Howitzer Component Mass Properties 18 

9. Length Properties 19 

B. EQUATIONS OF MOTION (EOM) 19 

1. Generalized Coordinates (q;) 22 

2. Kinetic Energy (T) 26 

3. Potential Energy (U) 28 

4. Generalized Power Lost Due to Damping Forces (F) 30 

5. External Forces (Qj) 31 

6. Assembling Lagrange's Equation 32 

C. HUMMV/HOWITZER NUMERICAL INTEGRATION 32 

1. Transformation to 1st Order Equations 32 

2. Integration Difficulties 34 

III. ROADWAY INPUTS 36 

A. WASHBOARD COURSE 36 

39 



B. "WORST CASE' POTHOLE COURSE 



IV. MODEL RESULTS 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOI 
MONTEREY CA 93943-5101 



.41 



A. WASHBOARD COURSE RESULTS 42 

1. Results at 3.75 MPH 43 

2. Results at 5.0 MPH 43 

3. Results at 5.25 MPH 43 

B. NATURAL FREQUENCIES AND MODE SHAPES 50 

C. "WORST CASE" POTHOLE COURSE RESULTS 50 



V. STRUCTURAL ANALYSIS 52 

A. A-FRAME CALIBRATION 52 

1. A-Frame Instrumentation 53 

2. Support Equipment 53 

3. Load Cells 56 

4. Calibration Procedures 56 

5. Calibration Factor 58 

6. A-Frame Stiffness 60 

7. Observations 61 

B. STRESS ANALYSIS 61 

1. T-Bar Stress 62 

2. A-Frame Stress 64 

3. Worst Case Stress and Factors of Safety 66 

VI. INSTRUMENTED TEST 67 

A. INSTRUMENTATION 67 

1. Accelerometers 67 

2. Single Strain Gages 69 

3. Other Gages 69 

B. TEST COURSES 72 

1. Belgian Block Course 72 

2. Six Inch Washboard Course 72 

3. Cross Country Course Number 1 73 

4. Cross Country Course Number 3 73 

C. TEST PROCEDURES 74 



v 



D. TEST RESULTS AND DISCUSSION 74 

1. Time History Results 75 

2. Frequency Domain 84 

3. Other Results 91 

VII. CONCLUSIONS/RECOMMENDATIONS 96 

A. CONCLUSIONS 96 

B. RECOMMENDATIONS 97 

APPENDIX A. SYSTEM PARAMETER VALUES 98 

APPENDIX B. MATRICES FOR EQUATIONS OF MOTION 101 

APPENDIX C. NUMERICAL INTEGRATION COMPUTER CODE 104 

APPENDIX D. POTHOLE COURSE COMPUTER CODE 108 

APPENDIX E. NATURAL FREQUENCIES AND MODE SHAPES 113 

APPENDIX F. A-FRAME CALIBRATION DATA 114 

APPENDIX G. HOWITZER ORIENTATION PICTURES AND 

DIAGRAMS 123 

LIST OF REFERENCES 127 



INITIAL DISTRIBUTION LIST, 



128 



LIST OF FIGURES 



2-1 HUMMV/Howitzer Model Components 6 

2-2 Force vs Velocity Plot for the Front HUMMV Shock 9 

2-3 Force vs Velocity Plot for the Rear HUMMV Shock 10 

2-4 Torsional Spring Stiffness for Serial No. 22 (Upward Motion).... 12 

2-5 Torsional Spring Stiffness for Serial No. 22 

(Downward Motion) 13 

2-6 Torsional Spring Stiffness for Serial No. 27 (Upward Motion).... 14 

2-7 Torsional Spring Stiffness for Serial No. 27 

(Downward Motion) 15 

2-8 Force vs Velocity Plot for the Howitzer Shock 16 

2-9 Length Properties for the HUMMV/Howitzer System 20 

2-10 Generalized Coordinates for the HUMMV/Howitzer System 23 

2- 11 Key Locations for the HUMMV/Howitzer System 25 

3- 1 Six Inch Washboard Course Wheel Inputs 38 

3- 2 "Worst Case" Pothole Profile and Axle Trace 40 

4- 1 Model Results for A-Frame Force at 3.75 MPH (Time History). .44 

4-2 Model Results for A-Frame Force at 3.75 MPH (FFT) 45 

4-3 Model Results for A-Frame Force at 5.0 MPH (Time History) 46 

4-4 Model Results for A-Frame Force at 5.0 MPH (FFT) 47 

4-5 Model Results for A-Frame Force at 5.25 MPH (Time History). .48 

4-6 Model Results for A-Frame Force at 5.25 MPH (FFT) 49 



Vll 



4- 7 Model Results for A-Frame Force for "Worst Case" Pothole 

at 15 MPH 51 

5- 1 Location of A-Frame Strain Gages 53 

5-2 Vertical Configuration for A-Frame Test 55 

5-3 Lateral Configuration for A-Frame Test 55 

5-4 A-Frame Connection to Test Mount 57 

5-5 A-Frame Apex Clamping Device 57 

5-6 Balance of Forces on the A-Frame 58 

5-7 T-Bar and A-Frame Assembly 62 

5-8 T-Bar Cantilevered Beam Model 63 

5- 9 Stress Concentration Due to a Circular Hole in a Pipe Section 65 

6- 1 Howitzer Instrumentation for APG Test 68 

6-2 Location of the Elevating Mass Accelerometer 70 

6-3 Location of T-Bar Strain Gages 71 

6-4 APG Belgian Block Course 72 

6-5 APG Six Inch Washboard Course 73 

6-6 A-Frame Force at 3.75 MPH on the Washboard Course 77 

6-7 Trail Vertical Acceleration at 3.75 MPH on the Washboard 

Course 78 

6-8 A-Frame Force at 5.0 MPH on the Washboard Course 79 

6-9 Trail Vertical Acceleration at 5.0 MPH on the Washboard 

Course 80 

6-10 A-Frame Force at 5.25 MPH on the Washboard Course 81 

viii 



6-11 Trail Vertical Acceleration at 5.25 MPH on the Washboard 

Course 82 

6-12 A-Frame Force at 15 MPHon Cross Country Course 

Number 3 83 

6-13 FFT of A-Frame Force at 3.75 MPH on the Washboard 

Course 85 

6-14 FFT of Trail Vertical Acceleration at 3.75 MPH on 

the Washboard Course 86 

6-15 FFT of A-Frame Force at 5.00 MPH on the Washboard 

Course 87 

6-16 FFT of Trail Vertical Acceleration at 5.00 MPH on 

the Washboard Course 88 

6-17 FFT of A-Frame Force at 5.25 MPH on the Washboard 

Course 89 

6-18 FFT of Trail Vertical Acceleration at 5.25 MPH on 

the Washboard Course 90 

6-19 Lateral Force from the APG Washboard Course at 3.75 MPH. ..92 

6-20 Lateral Force from the APG Washboard Course at 5.0 MPH 93 

6-21 Lateral Force from the APG Washboard Course at 5.25 MPH.. .94 

6-22 Lateral Force from Cross Country Course Number 3 at 

15 MPH 95 

G-l Towing in the Firing Position 123 

G-2 Towing in Stowed Position 124 

G-3 A-Frame/ T-Bar Support 124 

G-4 T-Bar (View from the Top Down) 125 

G-5 Howitzer Wheels and Suspension 126 



IX 



LIST OF TABLES 



A-l COMPUTER MODEL PARAMETER VALUES 99 

E-l NATURAL FREQUENCIES AND MODE SHAPES 113 

F-l VERTICAL TEST #1 A-FRAME CALIBRATION DATA 115 

F-2 VERTICAL TEST #2 A-FRAME CALIBRATION DATA 117 

F-3 LATERAL TEST #\ A-FRAME CALIBRATION DATA 119 

F-4 LATERAL TEST #2 A-FRAME CALIBRATION DATA 121 



x 



NOMENCLATURE 



SYMBOL 


DESCRIPTION 


HUMMV 


Abbreviation for the howitzer prime mover 


CG 


Center of gravity of component masses 


Ki 


The combined spring stiffness of the two front 
HUMMV tires 


k 2 


The combined spring stiffness of the two 
front HUMMV suspension springs 


k 3 


The combined spring stiffness of the two 
rear HUMMV tires 


K4 


The spring combined stiffness of the two rear 
HUMMV suspension springs 


K6 


The combined spring stiffness of the two 
howitzer tires 


k 7 


The combined spring stiffness of the 
two howitzer torsional springs 


K a 


The combined spring stiffness of the A- 
Frame and T-Bar 


Bi 


The combined damping constant for the two 
front HUMMV shock absorbers 


b 2 


The combined damping constant for the two 
rear HUMMV shock absorbers 


b 3 


The combined damping constant of the two 
howitzer shock absorbers 


B 4 


The rotational damping constant for the 
elevating clutch 



xi 



The combined mass of the two front HUMMV 
wheels and suspension components 

The combined mass of the two rear HUMMV 
wheels and suspension components 

The mass of the HUMMV body when it is 
fully combat loaded 

The mass of the howitzer wheel 

The total mass of the howitzer trail/saddle 

The mass of the howitzer elevating mass 

The mass moment of inertia of the HUMMV 
for pitch motion about its CG 

The mass moment of inertia of the howitzer 
trail/saddle mass about its CG 

The mass moment of inertia of the elevating 
mass about the trunion 

The static weight carried by the rear HUMMV 
wheels with the howitzer attached 

The static weight carried by the front 
HUMMV wheels with the howitzer attached 

The distance from the HUMMV CG 
to the connecting point of HUMMV front 
shocks and suspension springs 

The distance from the HUMMV CG 
to the connecting point of HUMMV rear 
shocks and suspension springs 

The distance from the HUMMV/howitzer hitch 
connection point to the HUMMV CG 

The distance between the center of the 
HUMMV front and rear axles 



The distance from the front HUMMV axle 
center to the howitzer axle center 

The distance from the howitzer axle to the CG 
of the howitzer trail/saddle mass (M5) 

The distance from the howitzer suspension 
spring and shock to the HUMMV/howitzer 
hitch connecting point 

The horizontal distance from the trunion to 
the A-Frame connection point on the T-Bar 

The horizontal distance from trunion to the 
HUMMV/howitzer hitch connection point 

The distance from where the point A-Frame 
attaches to the saddle, to the trunion 

The distance from CG of the gun tube to the 
trunion 

The static angle made between the horizontal 
plane and the gun tube 

The angle made between the A-Frame and 
the gun tube 

The vertical displacement for the CG of front 
HUMMV wheel mass 

The vertical displacement for the CG of rear 
HUMMV wheel mass 

The vertical displacement of HUMMV CG 

The vertical displacement of the CG of 
howitzer wheel mass 

The pitch angle of the HUMMV about its CG 



02 The pitch angle of the howitzer trail/saddle 

mass about its CG 

03 The rotation of gun tube about the trunion. 



xiv 



ACKNOWLEDGEMENT 



I thank Professor Healey for all his help in preparing this 
study. LTC Melvin Miller and Mr. Jerome Nathan of AMCCOM, Rock 
Island Arsenal, IL, the study sponsors, provided the funding to the 
Naval Postgraduate School for the study. CPT Dave Anderson's and 
CPT Steve Roger's efforts were instrumental in developing an original 
simplified model of just the howitzer, that served as the basis for the 
complete HUMMV/Howitzer system model I developed later and 
used to predict loads on the howitzer. 



I. INTRODUCTION 



A. GENERAL 

The purpose of this study was to determine the feasibility of 
transporting the Ml 19 howitzer in the firing position over terrain 
and at speeds used during normal tactical operations. The end 

objective was to determine if costly design changes would be 

required to safely permit towing the Ml 19 Howitzer in the firing 
position. The approach used was to model the vertical plane 

dynamics of the towing vehicle (HUMMV) and the howitzer, model 
representative terrain inputs to the wheels in the form of potholes, 

and to perform experiments to validate the model. It is proven that 
the simulation model is reasonable and it follows that stresses in the 
critical components are not expected to warrant redesign. 

B. BACKGROUND 

The Ml 19 Howitzer is a British designed howitzer that is 
currently being built at Rock Island Arsenal (RIA), Illinois and 
fielded to US Army Field Artillery units. The Ml 19 Howitzer 

operational manuals specify towing procedures that require the gun 

tube to be placed in the stowed position when moving the howitzer 

during tactical operations. These towing procedures are different 
and more difficult to implement than those used for all other towed 
howitzers in the US Army. All other US Army towed howitzers are 
towed during tactical operations with the gun tube in the firing 
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position, facilitating rapid emplacements and displacements that are 
required for successful tactical operations. 

The Ml 19 Howitzer has an existing A-Frame brace which can 
support the gun tube when the howitzer is towed in the firing 
position. The primary concern of this study was whether or not the 
A-Frame and the T-Bar structure on the howitzer’s trail section, to 
which the bottom of the A-Frame brace connects, will be strong 
enough to support the loads generated when towing the howitzer 
over rough terrain. This study also addresses concerns about trail 
loads on the howitzer. This study does not evaluate rollover stability 
for the Ml 19 Howitzer when towing in the firing position. 

The terrain that the HUMMV/Howitzer system can traverse is 
clearly limited by the terrain the HUMMV can negotiate, therefore, 
the Ml 19 Howitzer and HUMMV have been evaluated as a system. A 
dynamic simulation model of this system was seen to be a necessary 
tool to aid in making design decisions from a rational viewpoint and 

is the basis of this study. 

C. METHOD OF ANALYSIS 

To analyze the Ml 19 Howitzer a computer model of the 

HUMMV/Howitzer system was developed to predict the response of 
the system to roadway inputs. This computer model treats the major 

mass elements of the HUMMV/Howitzer system as lumped masses to 
develop vertical plane heave and pitch equations of motion for the 
system. The development of this computer model is contained in 
Chapter II. To drive the computer simulation model, terrain models 
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were developed for both a "worst case" trench type pothole and a 
sinusoidal "washboard" roadway. Details on the development of 

these roadway inputs are contained in Chapter III. Chapter IV 
contains the results of driving the computer model over the roadway 
inputs. Chapter V explains the structural analysis that was used to 
determine the maximum stress in the A-Frame and the T-Bar. 

Chapter V also contains the procedures and results of calibrating an 
A-Frame that was used as a load cell for an instrumented dynamic 
test. Chapter VI explains the conduct and the results of the 
instrumented test on the Ml 19 Howitzer/HUMMV system that was 
conducted at Aberdeen Proving Grounds, Maryland in April 1993. 

Also, in Chapter VI the results of the model are compared to the 
experimental results from the instrumented test. This comparison 
showed reasonable correlation between the model and experimental 
results. Conclusions and recommendations based on this study are 
contained Chapter VII. The final recommendation is that towing in 
the Firing position should be allowed. 
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II. IIUMMV /HOWITZER SYSTEM COMPUTER MODEL 



The computer model developed is based only on the heave and 
pitch equations of motion (EOM) for the HUMMV/Howitzer system. 
The reason for this simplification is that the lateral, yaw, and roll 
motions are much less significant for determining the forces in the 
A-Frame for the type terrain inputs considered by this study. While 
horizontal plane motions are key elements for the analysis of turns 
and rolling during transit, the concerns of this study indicate that a 
vertical plane model is adequate. The major mass groups of the 
HUMMV/Howitzer were treated as lumped masses in the 
development of the equations of motion. The equations of motion 
were developed using Lagrange techniques. 

A. HUMMV/HOWITZER SYSTEM PARAMETERS 

The primary sources of information for system parameters for 
the HUMMV and howitzer were the US Army Tank-Automotive 
Command Research, Development & Engineering Center Technical 
Report No. 13337, [Ref. 1], for HUMMV data and US Army Combat 
Systems Test Activity, APG, Report No. P-83179, [Ref. 2], for the 
Ml 19 Howitzer data. Additionally, many system parameters were 
determined from direct measurement on the HUMMV and howitzer 
or from design diagrams, or obtained from the study sponsor at Rock 
Island Arsenal (RIA). For readers of this study that are unfamiliar 
with the components of the Ml 19 Howitzer, Appendix G contains 



4 



several diagrams that depict the howitzer components that are 
important for this study. 

The system parameters used for the HUMMV are based on using 
the Model No. Ml 069 HUMMV as the prime mover for the Ml 19 
Howitzer. The cargo loading of the HUMMV significantly effects its 
center of gravity and consequently effects all of the HUMMV's 
system parameters which are dependent on its center of gravity. 

Modeling was therefore limited to using the same cargo loading 
configuration that was used in [Ref. 1] and other loading effects were 
ignored. These HUMMV simplifications are justified because for this 
computer model the HUMMV serves primarily as the driving force at 
the howitzer hitch causing pitching motions of the entire howitzer 
and we are not concerned with the loading of HUMMV components. 

The major components of the HUMMV/Howitzer system computer 
model are identified in Figure 2-1 below. 

In the model, the tires act as linear springs which provide force 

inputs to the HUMMV and howitzer wheels, based on the terrain over 
which the model was "driven". The suspension systems for both the 
HUMMV and the howitzer are modeled as springs and dampers 
connected in parallel. The A-Frame and T-Bar are modeled as a set 

of stiff springs connected in series. The following paragraphs 
provide more detail on the modeling of the individual components 

and the determination of appropriate parameter values. Appendix A 

provides a listing of the numerical values, units, and the symbols 
used for each system parameter in the computer model. The 
parameter values for those components with identical components on 
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Figure 2-1 HUMMV/Howitzer Model Components 
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the right and left sides of the system (wheels and suspension 
components) are multiplied by two. The values in Appendix A for 
these components reflect the effects of both sides. 

1. HUMMV and Howitzer Tires 

The HUMMV has the same tires on the front and rear, 

however, the front and rear tire pressure is different (20 psi vs 30 
psi) resulting in the rear tire spring being much stiffer. The 
stiffness used for the HUMMV tires was obtained from normal force 
vs tire deflection data [Ref. 1]. The stiffness of the howitzer tire used 
in the model was estimated because no measured data on the 
howitzer tire stiffness was available. 

The tires are modeled as linear springs which transmit a 

vertical force to the respective wheel masses proportional to the 
amount that the springs are stretched or compressed. The amount of 
force is determined by the relative dynamic displacement between 
the terrain input for the wheel and the center of mass of the wheel. 
If the tire leaves contact with the ground there will be no input from 
the terrain until it again returns to contact. To account for this loss 
of contact situation in the computer model, the tire spring force is set 
to the static weight that the tire supports when the tensile (de- 
compressive) spring force exceeds the static weight on the tire. The 

tire springs are turned back on in the computer model when the 
conditions are met for road surface contact to be re-established. 

2. HUMMV Suspension System 

The HUMMV suspension system consists of a shock 
absorber and coil spring for each wheel. The shock absorber passes 
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through the middle of the coil spring. The force in the spring is 
linearly proportional to the relative displacement between the center 

of mass of the wheel and the attachment point on the HUMMV body. 
The values for the spring constants were obtained from Table 5-13 
of [Ref. 1]. Precise modeling of the shock absorbers is difficult 

because they have a mid-stroke region and a hydraulic bump stop 

region. In the mid-stroke region the shock absorber force is linearly 
dependent on the relative velocity of the center of mass of the wheel 
and the attachment point on the HUMMV. However, the linear 

damping constant changes depending on whether the force in the 

shock absorber is tensile or compressive. When the shock absorber 

is compressed or stretched into the hydraulic bump stop region the 
damping of the shock absorber greatly increases. This sudden 

increase in damping causes problems in the integration routine that 
was used for the computer model. To overcome these problems and 

still account for the energy loss due to the increased damping in the 
hydraulic bump stop region, the mid-stroke linear damping constant 

obtained from [Ref. 1] was increased by 20% for both the front and 

rear shock HUMMV absorbers. Figures 2-2 and 2-3 shows the force 
versus velocity plots that were obtained from [Ref. 1] for the front 

and rear shock absorbers, respectively. 

3. Howitzer Suspension System 

Modeling the howitzer suspension system is more 

complicated than the HUMMV because it consists of a torsional spring 
and a shock absorber that is not vertically oriented. The torsional 
spring is modeled as a linear vertical spring that is connected to the 
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Figure 2-2 Force vs Velocity Plot for the Front HUMMV Shock 
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Figure 2-3 Force vs Velocity Plot for the Rear HUMMV Shock 



center of mass of the howitzer wheel and attached to the howitzer 
trail directly above the wheel center of mass. The effective spring 

stiffness constant was calculated based on the measured data of the 

torsional spring stiffness. Figures 2-4 through 2-7 show the angular 

deflection versus vertical force plots that were contained in [Ref. 2]. 
Using the average torsional spring stiffness from this data and the 
length of the road arm, an angular rotation was converted to an 
equivalent vertical displacement at the end of the road arm. The 

spring constant for the model was then calculated by determining 
the constant value that must be multiplied times this vertical 

displacement at the end of the road arm, to produce a moment that is 
equivalent to that caused by the torsional spring for a specified 

angular rotation. 

A similar procedure was used to model the howitzer's 

shock absorber as a vertical shock absorber with a linear damping 
constant that is velocity dependent. The attachment points for the 

ends of the shock absorbers in the model are the same as the 

modeled end points of the howitzer suspension spring. The actual 

non-vertical shock absorber force was transformed to a vertical force 
in the model by determining the amount of vertical force applied a t 
the end of the road arm needed to produce an equivalent moment 

about the road arm pivot point. Like the HUMMV shock absorbers, 
the howitzer shock absorbers are non-linear and have different 

damping values for tension and compression that are accounted for 
in the computer model. This non-linear damping in the shock 

absorbers is a significant aspect of the computer model. The values 



used for the shock absorber damping constants were determined 
from force versus velocity test data provided by RIA. Figure 2-8 
shows the force versus velocity plot for the howitzer shock absorber. 




co 



o 

CO 



co 



o 

CO 



on 

<U 

u 

-a 

I 



o 



a 

0) 

_ C 

uo cj 

^ a 

u 

jjj 

a 

taO 
O 2 
— < 






o 



Figure 2-4 Torsional Spring Stiffness for Serial No. 22 (Upward 
Motion) 



Ml 19 HOWITZER, 105mm GUN, Vehicle No. 22 

Date of Test: 4 through 5 April 1989 




Figure 2-5 Torsional Spring Stiffness for Serial No. 22 (Downward 
Motion) 
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Figure 2-6 Torsional Spring Stiffness for Serial No. 27 (Upward 
Motion) 
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Figure 2-7 Torsional Spring Stiffness for Serial No. 27 (Downward 
Motion) 
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Figure 2-8 Force versus Velocity Plot for Howitzer Shock 



4. A-Frame/T-Bar Stiffness 

The A-Frame and T-Bar are modeled as stiff springs 
connected in series to determine a combined linear spring stiffness. 
The stiffness of the A-Frame was determined during the A-Frame 
calibration test (explained in Chapter V). The stiffness of the T-Bar 
was determined by using the principle of superposition of forces to 
determine the effective bending stiffness for the non-symmetric 
bending of the T-Bar that is caused by the A-Frame. This non- 
symmetric bending of the T-Bar results because the angle at which 
the A-Frame transmits forces and moments to the T-Bar is not 
perpendicular to an axis of symmetry of the T-Bar. 

5. Equilibrator Springs 

The equilibrator springs provide only small dynamic forces 
because their combined stiffness is nearly 1000 times less than the 
effective stiffness of the A-Frame/T-Bar assembly. Therefore, 
equilibrator springs are neglected in the computer model. 

6. Elevating Clutch 

The effects of the elevating clutch of the howitzer were 
neglected in the final computer model. It is believed backlash in the 
elevating mechanism gears will prevent the clutch from being 
activated while towing with the A-Frame attached because the 

relative rotation of the gun tube is very small. However, because of 

some uncertainty about the clutch several computer runs were made 
with clutch modeled as a linear torsional damper. The value of the 
torsional damper used in the model was estimated because no test 
data on the clutch was available. These model runs with the clutch 



included, produced only slightly different results from the runs 

without the clutch included. Therefore, because of the previously 
mentioned belief that the clutch would not be activated, it was 
neglected in the final model. 

7. HUMMV Component Mass Properties 

The HUMMV is divided into three lumped masses, the 
HUMMV body mass and two wheel masses. The front and back 

HUMMV wheels are treated as having identical mass values in the 
model. Each wheel weight/mass includes the weight of the tire, 
wheel, shocks, suspension springs, and suspension support arms. The 

HUMMV body weight/mass varies dependent upon it's cargo load. 
The location of the center of gravity also can vary and is determined 
by the positioning of the cargo. Therefore, the values used in the 

computer model for the HUMMV mass, the mass moment of inertia 
for pitch, and lengths that are based dependent on the center of 
gravity, are estimated values. These estimated values are based on 

the same loading situation as was used in [Ref. 1]. The pitch mass 

moment of inertia for the HUMMV is about the center of gravity of 

the sprung HUMMV mass. 

8. Howitzer Component Mass Properties 

The howitzer is divided into three separate lumped masses 

in the model; the howitzer wheel masses, the elevating mass, and 

the combined trail/saddle mass. The wheel weight/mass includes 
the tire, rim, and approximately 30% of the road arm that rests on 
the wheel assembly. The elevating mass includes the gun tube, the 
cradle, the recoil system and all other parts that move the when the 



gun is elevated about the trunion. The mass moment of inertia of the 
gun tube for rotation about the trunion was estimated using the 
approximate masses of all the elevating components. The center of 
gravity location for the elevating mass was estimated by RIA. The 
combined trail/saddle mass includes all of the gun mass except the 
elevating mass and the wheel mass. The center of gravity location 
for the trail/saddle mass and the pitching mass moment of inertia for 
the trail/saddle mass was estimated from design diagrams. 

9. Length Properties 

All lengths used in the computer model are identified in 

Figure 2-9 below. These length properties were determined either 

by direct measurements on the howitzer or by scaling from design 
diagrams. 

B. EQUATIONS OF MOTION (EOM) 

In this section the heave and pitch equations of motion (EOM) 
for the model are developed from Lagrange's equations using the 
HUMMV/Howitzer system parameters that were previously 
explained. Lagrange’s equations require the determination of the 
system’s total kinetic and potential energy and the power dissipated 
by the dampers, in terms of a set of generalized coordinates. 

Throughout this development of the EOM, expressions are 

developed that are functions of generalized positions and velocities, 
angular displacements, angular velocities, and angular accelerations. 
To simplify the EOM with minimal loss of representational accuracy, 

these expressions are linearized using small angle approximations. 
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2-9 Length Properties of HUMMV/Howitzer System 
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The sine of a small angle is approximated as just the angle (sin0=0) 
and the cosine of a small angle is approximated as equal to one 
(cos0=l). These approximations are valid because the pitch or 
angular rotations are predicted to be less than five degrees. Order of 

magnitude considerations are used to further simplify the EOM. 
Because the angular displacements are small their angular velocities 
and angular accelerations are also small. Therefore, any terms that 
contain an angular displacement, angular velocity, or angular 
acceleration multiplied times another angular displacement, angular 
velocity, or angular acceleration are neglected because these terms 
will be very small. 

The general form of Lagrange's equations used is: 



qi = generalized coordinates 
T = total kinetic energy of the system 

U = change in the potential energy of the system with respect to 
its potential energy in the static-equilibrium position 
F = the power dissipated by forces in dampers 
Q = the input force acting on the system. 

The development of the above components of Lagrange’s equation 

and the assembly of these components into the system’s heave and 

pitch EOM follows. Due to the length of this development the reader 
is provided only the final results from each step. The intermediate 

multiplication's, differentiation's, and linearization's are omitted. 



d_r dT 

dtU<ii 



dT dU dF 

dq ; dq . dq. 



where: 



1. Generalized Coordinates (qj) 

The set of generalized coordinates, q ; 5 that was used to 
.uniquely specify, the motion of the HUMMV/Howitzer system are 

depicted on Figure 2-10 and are defined follows: 

* Z 1 &Z 2 ' The vertical displacements for the front and rear 

HUMMV wheel masses, respectively. 

* Z 3 : The vertical displacement of the center of gravity of the 

HUMMV body. 

* Z 4 : The vertical displacements of the center of gravity of the 

howitzer wheel mass. 

* 0i.' The pitch angle of the HUMMV body about its center of 

gravity. 

* 63 : The pitch angle of the howitzer trail/saddle mass about its 

center of gravity. 

* 63 : The absolute rotation of the gun tube. 

As depicted on Figure 2-10, upward vertical motions and clockwise 
rotational motions are defined as positive. The equilibrium position 

for each of the generalized coordinates is the static equilibrium 

position for the respective components that are correlated to each of 
the generalized coordinates. In defining positions to key locations in 
the system unit vector notation of i, j, and k will be used 

corresponding to the x, y, and z axis directions of the Cartesian 

coordinate reference system shown on Figure 2-10. 

The kinetic and potential energy of the system and the 
damping and input forces are all computed in terms of these 
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Figure 2-10 Generalized Coordinates for the HUMMV/Howitzer 
System 



DEFINITION OF WAY INPUTS- 
u i : The road input for the front HUMMV wheel 
U2: The road input for the rear HUMMV wheel 
U3 : The road input for the howitzer wheel 



generalized coordinates and the road inputs Ui, U 2 , and U 3 . To carry 
out the required computations, position vectors to key locations o n 

the HUMMV and howitzer are defined in terms of the generalized 
coordinates. The position vectors to these key locations are depicted 
on Figure 2-11 and defined as follows: 

* rj & ti = the front and rear attachment points, respectively, for 

the suspension springs and shock absorbers on the HUMMV. 

* r 3 = the center of mass of the trail/saddle mass 

* r 4 = the modeled attachment location of the suspension spring 
and shock absorber of the howitzer 

* 1-5 = the location where the A-Frame attaches to the T-Bar 

* i *6 = the location of the center of mass of the gun tube/elevating 

mass 

* r 7 = the location where the A-Frame attaches to the gun 
tube/elevating mass 

These key locations that are defined above can be expressed as 

position vectors in terms of the above defined lengths (Figure 2-9), 
generalized coordinates (Figure 2-10), and the unit vectors i, j, and k 
of the Cartesian coordinate reference system (Figure 2-10) as follows: 

* r, = (Z 3 - LiSin(0,))k + (L 1 cos(0 1 ))i 

* r 2 -(Zj + L 2 sin(0 1 ))k - (L 2 cos(0j))i 

* r 3 = (Z, + LjSin^,) + L u sin(0 2 ))k - (LjCos^,) + L u cos(0 2 ))i 

* r 4 = (Zj + LgSin^j) + L 12 sin(0 2 ))k - (L 3 cos(0 1 ) + L, 2 cos(0 2 ))i 

* r 5 =(Z 3 + LjSin^) + (1^+ L 16 )sin(0 2 ))k - (L 3 cos(0,) + (L^ +L 16 )cos(0 2 ))i 
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Figure 2*11 Key Locations for the HUMMV/Howitzer System 
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= the location where the A-Frame attaches to the T-Bar 
= the location of the center of mass of the gun tube/elevating mass 
= the location where the A-Frame attaches to the gun tube/elevating mass 



* 



* 



= (Zj + Lj sin^j) + L 16 sin(0 2 ) + L CG sin(0 o + 0 3 ))k - 
= (Z 3 + L 3 sin(0 1 )+L 16 sin(0 2 ) + L A sin(0 o +0 3 ))k - 



r L 3 cos^) + L 16 cos(0 2 )Y 

+ ®j) ) 

^L 3 cos(0,) + L 16 cos(0 2 )Y 
^+L a cos(0„ + 0 } ) / 



The velocities of these key locations can be determined by 
taking the time derivatives of the position vectors. However, only 
the velocities of locations where there are velocity dependent actions 
occurring in Lagrange's equations need to be determined. The 
required velocities of the key locations on the HUMMV and howitzer 
are as follows: 



* r, = (z 3 - L 1 cos(0 1 )0,)k - (Lj sinCO^O ,)i 

* r 2 = (Z3 + L 2 cos(0,)0,)k + (L 2 sin(0 1 )0,)i 

* r 3 = (z, + Ljcosf©,)©! + L n cos(0 2 )0 2 )k + (LjSin^)©! + L u sin(0 2 )0 2 )i 

* r 4 = (Zj + Ljcost©,)©] + L 12 cos(0 2 )0 2 )k + (L 3 sin(0 1 )0, + L 12 sin(0 2 )0 2 )i 



r 6 = 



Z 3 + L 3 cos(0,)0, + L 16 cos(0 2 )0 2 
+ Lpq cos(0 o + 0 3 )0 3 



k + 



L 3 sin(0 1 )0 1 + Lj 6 sin(0 2 )0 2 
k. ^cg sin (0 o + 0 3 )0 3 



In addition to the above velocities, the vertical velocity of the CG of 
the HUMMV, Z 3 , and each of the wheel masses, Z p Z 2 , & Z 4 are 

needed in Lagrange's equation. 

2. Kinetic Energy (T) 

The velocities defined above are used to calculate the 
kinetic energy of the system. The kinetic energy is calculated by 
summing the linear and rotational kinetic energy for each of the 



system's masses. The total kinetic energy of the system is given by: 
I T, - £ (*m,v, 2 ♦ 1*1,6/) i - 1,6; j - 1,3 

where: 

T; = the kinetic energy of each mass in the system 
m j = the mass of each lumped mass 
Vi = the linear velocity of the each mass 

I j = the heave mass moment of inertia of each rotating mass 
about its center of gravity 

0j = the angular velocity of each rotating mass 

Therefore, the total kinetic energy of the HUMMV/howitzer system 

was determined by summing the kinetic energy components as 

follows: 

T=^M 1 (Z 1 ) 2 + 4M 2 (Z 2 ) 2 + jM i (Z 3 ) 2 + tM 4 (Z 4 ) 2 + 4M 5 (r 3 ) 2 

iMdO 2 + iW ) 2 + tM ^) 2 + ii ,(«,) 2 

In the above expression the square of the velocities is determined by 

substituting in the appropriate velocity expression as defined above, 

and taking the dot product of the velocity with itself. This 

substitution and dot product operation is left for the reader. When 

taking this dot product all of the k dot i and i dot k cross product 

terms drop out because the i and k unit vectors are orthogonal. 

Having determined the total kinetic energy, the time derivative of 

the partial derivative of the total kinetic energy with respect to the 

d dT 

velocity of each of the generalized coordinates, —77 can be 

dt dg t 

determined. Carrying out this differentiation and linearizing yields 
the following: 



* !£L 

dt dZ, 
. d dT 
dt dZ 2 
d dT 

* dt dZ 3 



M,Z, 

M 2 Z 2 

(M 3 + M s + M 6 )Z 3 + (M 5 L 3 + M 6 L 3 )0 1 + (M 5 L u + M 6 L 16 )0 2 



+ M 6 L CG cos(0 o )0 3 



, d dT 
dt dZ 4 
d dT 
* dt d0, 



m 4 z 4 

(M s L 3 + M 6 L 3 )Z 3 + (I, + M s L 3 2 + 



M 6 L 3 2 )0, 



+ (M s L 3 L n + M 6 L 3 L 16 )0 2 + MgLjLpQ cos(0 o )0 3 
(M. 1 -.. + M S L .^Z, + (M 5 L 3 L u + MjLjL.je, 

* dt ou 2 

+ (I 2 + M 5 L h 2 + M 6 L 16 2 )0 2 + M 6 L 16 L CG cos(0 o )0 3 

d dT 

M 6 L cg cos(0 o )Z 3 + M 6 L 3 L CG 005(00)0, + M 6 L, 6 L cg cos(0 0 )0 2 

* dt 50 3 

+ (I 3 + M 6 L cg )0 3 



The total kinetic energy, T, must also be differentiated 

with respect to each of the generalized positional coordinates, 

however, for this system this differentiation is equal to zero for all 

coordinates. 

3. Potential Energy (U) 

The potential energy for this system results from the 
strain energy that is developed in the suspension springs and the A- 

Frame spring. The gravitational potential energy that arises from 

changes in vertical positions is ignored because the resulting terms 
drop out of the dynamic equations of motion. The potential energy 
in any spring is determined as follows: 



2 8 



Uj = 2^i(AXj) 2 



where: 

U; = the potential strain energy in each spring 

K; = the spring stiffness for spring 

AX; = the change in length of the spring. 



Since locations to the attachment points of the ends of all the springs 
have been defined above in a common inertial coordinate reference 
system, taking the difference of the position vectors of the end points 
of each spring will yield the change in length of the spring. Thus for 

this system, the total potential energy is found by summing the 
potential energy in each of the springs as follows: 




* I K 4( r 2 - 
* T^a( f 7 " r s) 




Therefore, the partial derivative of the total potential energy with 

respect to each of the generalized coordinates, after linearizing 

dq ; 



results in the following: 



* au 

* 

dZ, 

, eu 
az 2 
eu 
* az^ 



= k,z, + k 2 z, - k 2 z 3 + k 2 l, 0, - K^Uj 

= K 3 Z 2 + K 4 Z 2 - K 4 Z 3 - K 4 L 2 0! - k 3 u 2 
= -K 2 Z, - K 4 Z 2 + (K 2 + K 4 + K 7 )Z 3 - k 7 z 4 



+ (K 2 Lj + K 4 L 2 + K 7 L 3 )6j + K 7 L, 2 0 2 

* = -K 7 z 3 + (K 6 + K 7 )Z 4 - K.L 30 J ~ K 7 L 12 0 2 - K 6 U 3 

o Z 4 
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k 4 l 2 z 2 + (-k 2 l, + k 4 l 2 + K 7 L 3 )Z 3 - k 7 l, 



au 
* as, 



k 2 l,z, - 



(K 2 L, 2 + K 4 L 2 2 + K 7 L 3 2 )S, + K 7 L 3 L 12 6 2 
au = K 7 L 12 Z 3 - K 7 L 12 Z 4 + k 7 l 3 l,^, + K 7 L 12 2 e 2 



* as, 



+ K a L a L 13 cos(S 0 )S 2 K a L a L 13 cos(S 0 )S 3 

au 

* — K A L A L 13 cos(S 0 )S 2 + K a L a L, 3 cos(0 o )0 3 

o o 3 



4. Generalized Power Lost Due to Damping Forces (F) 

The amount of generalized power that is dissipated by 
each of the dampers in the system given by the equation: 



F, = iB .(AkJ 



where: 

Fj = the generalized power dissipated in each damper 

Bj = the damping constant for each damper 

Ax, = the difference in velocity of the attaching points of each 

damper. 



The total amount of generalized power dissipated in the system is 
found by summing the energy dissipated in each dampers. Summing 

the energy dissipated by the damping forces yields the following: 



F=iB,(r, - Z J + - Z 2 ) ! ♦ - zj * - 6 2 f 



In Lagrange's equation this total dissipated energy, F, must be 

differentiated with respect to the velocities of each of the generalized 
dF 

coordinates, ■— . 

aq, 
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Carrying out this differentiation and again linearizing results in the 
following: 



dZj ^1^3 + BjLjOj 

aF 

B 2 Z 2 B 2 Z 3 B 2 L 2 0j 
aF • 

0 %-- -B 1 Z 1 - B 2 Z 2 + (B 1 +B 2 + ^)Z 3 - B 3 Z 4 + (~ B, L 3 + B 2 L 2 )0, 

aF 

^3^4 BjZj ~ ^3^3®1 ~ ^3^12^2 



4 



aF 
* ae, 



-«q — B 2 L 2 Z 2 + + B 2 L 2 + B 3 L 3 )Z 3 B 3 Lj 



+ (B.L , 2 + B 2 L 2 2 + B 3 L 3 2 )e; + B 3 L 3 Lj ^ 2 

aF • , 

B 3 Lj 2 Z 3 - B 3 L 12 Z 4 + BjLjLj^j + (B,L 12 + B 4 )0 2 — B 4 0 3 

* W } = ~ B ^ 2 + B 4 ® 3 



5. External Forces (Qi) 

The external forces acting on the system are the roadway 
inputs that are transmitted through the tire springs. For the 

computer model these input forces are treated as purely vertical 
forces. The forces result from the road inputs causing a change in 

length of the tire springs. These road inputs are depicted on Figure 
2-10 and defined as follows : 

* Ui = the road input for the front HUMMV wheel 

* U 2 = the road input for the rear HUMMV wheel 

* U 3 = the road input for the howitzer wheel. 



6. Assembling Lagrange's Equation 

Combining the parts of Lagrange's equations that were 
developed in the above paragraphs results in a system of equations 
of motion in the form: 

[M]|x| ♦ [C(x)]{x| ♦ [K]|x| - [f] 

where: 

[M] = the 7x7 system mass matrix 

[C(x)] = the 7x7 system damping matrix that is a function of 
whether the shock absorbers are in tension or compression 
[K] = the 7x7 system stiffness matrix 

{x}, {x}, & {x} = 7 x 1 column vectors that contain accelerations, 
velocities and positions of all coordinates, respectively. 

[f] = the 7x1 forcing vector that results from the road inputs. 

The [M], [C(x)], and [K] matrices for this system are contained in 
Appendix B. This completes the development of the heave and 

pitch EOMs. The resulting EOMs were numerically integrated to 

determine the time history for motions of the component masses of 
the system. 

C. HUMM V/HOWITZER MODEL NUMERICAL INTEGRATION 
The following section describes the procedures that were used 
to numerically integrate the previously developed equations of 
motion. 

1. Transformation to 1st Order Equations 

In order to numerically solve the equations of motion, the 
formulated system of second order ordinary differential equations 
must be transformed into a system of first order differential 



equations. This equation transformation was performed using matrix 
algebra and results in the following expression for the system of first 
order differential equations: 



where: 



[xdot] 



r [o] m i 


"x 




'0 ' 


Lf-M-'K] [~M~ l C(x)J J 


X 


+ 


_M _1 f_ 



[xdot] = a 14 x 1 column vector of which the first seven elements 
are the velocities of all the generalized coordinates and the last 
seven elements are the accelerations of the generalized coordinates. 



x 

x 



= a 14 xl column vector of which the first seven elements 



are the positions in terms of the generalized coordinates and the last 
seven elements are the velocities of the generalized coordinates. 



[0] = a 7 x 7 zero matrix 



[I] = a 7 x 7 identity matrix 

[-M _1 C(x)] & [-M’‘k] = these are 7 x 7 matrices that result from 

multiplying the negative inverse of the system mass matrix, [M], 
times the damping matrix, [C], and the stiffness matrix, [K], 
respectively 



0 

M''f 



= a 14 x 1 vector which the first seven elements are zero and 



the last seven elements are the matrix product of multiplying the inverse 
of the mass matrix, [M], times the forcing vector, [f|. 



For the integration the system was started from rest, thus, all 

seven of the system's displacement, {x}, and velocity, {x}, state variables 

had a value of zero for initial conditions. 



2. Integration Difficulties 

The equations of motion were numerically integrated using the 
Matlab ODE 45 integration routine. This is a variable time step 4th order 

Runge-Kutta integration routine. Integrating these equations of motion 

using ODE 45 was complicated by the non-linear damping of the shock 
absorbers and the required shutoff of the tire springs when they leave 
contact with the road surface. 

When the damping values of the shock absorbers would abruptly 
change from compression to tension and vice versa, there would be a 
corresponding abrupt change in the shock absorber force. Determining 

when the shock absorber damping value should change is built into the 
integration code by checking for each time step if each of the shock 

absorber forces are compressive or tensile. An "if" statement ensures that 
the correct damping value is used for the condition that exists. The abrupt 
changes in force in the shock absorbers slowed down the tolerance based 

ODE 45 routine. A similar situation occurred when the tire spring forces 
would abruptly change to the respective static weights when contact with 
the road was lost for a particular wheel. 

To increase the speed of the ODE 45 routine the integration was 
divided into one second time segments so the solution matrix that contains 

the state variables remains relatively small. The displacements and 

velocity states at the end of the previous time segment are used as the 
initial conditions for the next one second time segment. 

The resulting time histories for the displacements and velocities 
allowed post-processing to be completed determine accelerations, forces, 



and stress in key system components. The computer code developed for 
this numerical integration is contained in Appendix C. 



III. ROADWAY INPUTS 



This chapter explains the road inputs that were used to drive 
the computer model. The roadway inputs are defined in terms of the 
vertical change upward or downward of the values of the road inputs 
Ui, U2, and U3 that are shown on Figure 2 - 10 . These are independent 
vertical displacement inputs to tires for each of the wheels in the 

computer model. As defined in Chapter II, there is one input for the 
front HUMMV wheels, (Ui), a second input for the rear HUMMV 
wheels, (U2), and a third input for the howitzer wheels, (U3). The 
model can be driven over any terrain that can be defined in terms of 
the road inputs U2, and U3. The road inputs are functions of time 
which allows the user to control the speed at which the model is 
driven by using the time, velocity, and distance relationships. For 
this design study the road inputs used were a washboard course that 
was like one of the courses used in the actual test at Aberdeen 
Proving Grounds, MD (see Chapter VI), and a predicted "worst case" 
pothole road. 

A. WASHBOARD COURSE 

The washboard course that was used to drive the model was 

constructed to be similar to the Six Inch Washboard Course at APG. 

During the instrumented test at APG there were several data runs 
taken on the Six Inch Washboard Course (again, see Chapter VI). 
Comparing the model results for this type of washboard course to the 
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actual instrumented data from the APG test for a similar washboard 
provided a means of validating the output results of the computer 

model. 

To model the APG Six Inch Washboard Course a sinusoidal input 
was used with an amplitude of three inches and a distance from peak 
to peak of 72 inches. However, because the distance between the 
front HUMMV wheels and both the rear HUMMV wheels and also the 
howitzer wheels does not correspond to even multiples of the 72 inch 
peak to peak distance, a lag had to be built into the input sinusoid for 
both the rear HUMMV wheels and the howitzer wheels. This lag was 
necessary, to get the correct timing of the bump wheel inputs for all 
the wheels in the system. 

Additionally, in order to prevent difficulties in the integration 
routine that would occur if the wheel inputs were all turned on at 
the start of the integration run, the rear HUMMV and howitzer wheel 
inputs were not turned on initially. These inputs were turned on at 
the time when the respective wheels would start hitting the 
washboard course bumps if the HUMMV/howitzer system entered 

the course from level terrain. In other words, the front HUMMV 

wheels entered the washboard and started taking road input, Ui, 
while the rear HUMMV wheels and the howitzer wheels remained on 
level terrain which provided zero input to them. Then the rear 

HUMMV wheels entered the washboard and started taking road 
input, U 2 , while the front HUMMV wheels continued to have the 
input, Ui, and the howitzer wheels continued to have zero input. 



Finally, the howitzer wheels entered the washboard and at this point 

all three wheels had their separate inputs Ui, U 2 , and U 3 . 

This sinusoidal input, for which the wavelength was defined in 
terms of horizontal distance between peaks, was converted into a 
function of time by using time, velocity, and distance relationships. 
To calculate the time between successive peaks for a specified 
velocity the inputs Ui, U 2 , and U 3 were then defined in terms of time 
which allowed the computer model to be "driven" at any desired 

speed. Making the input a function of time was required for using 

the time based ODE 45 integration routine of Matlab to solve the 

equations of motion. Figure 3-1 depicts the wheel inputs U 2 , and 
U 3 for the washboard course being driven at 5.0 MPH. 



Washboard Wheel Inputs 
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Figure 3-1 Six Inch Washboard Course Wheel Inputs 



B. "WORST CASE" POTHOLE COURSE 

The model was also "driven" over a trench type pothole to 

simulate some of the most severe road conditions that the howitzer 
would be subjected to during tactical operations. The model results 
for A-Frame force from this terrain were used to predict if the A- 

Frame and T-Bar were strong enough to handle this "worst case" 

terrain. The pothole course road inputs, Ui, U 2 , and U 3 , are 
determined by a separate computer program (see Appendix D) that 
computes the changes in axle height as each of the HUMMV wheels 
and the howitzer wheels roll through the pothole. The road inputs 

are based on the path that would be traced out by the center of a 17 
inch radius wheel rolling through the specified pothole. The 

computer program that was developed to determine the road inputs 

allows the user to specify the desired dimension of the pothole for 

driving the computer model. For this design study the "worst case" 
pothole was determined to be 20 inches deep and 30 inches across 
with 85 degree slopes on each edge. This size was selected because it 
was judged to be the largest pothole that the HUMMV driver would 
traverse without stopping the HUMMV. Thus, the "worst case" 
pothole dimensions were selected based on what was judged to be 

the most severe terrain the HUMMV/howitzer would encounter 
while traveling during tactical operations. The "worst case" pothole 
selected was slightly larger than the dimensions of the pothole that 
was used in an informal test conducted at RIA in June 1992. Figure 
3-2 shows the pothole profile, and the axle path traced for the 
howitzer or HUMMV wheel rolling through this "worst case" situation. 




Figure 3-2 



"Worst Case" Pothole Profile and Axle Trace 



IV. MODEL RESULTS 



This chapter contains the results of driving the computer model 
over the terrain that was explained in Chapter III. The results from 

driving the model over the washboard course show close agreement 
actual results obtained from the APG test. A detailed comparison of 
the model results to the actual results is contained in Chapter VI. 

The model's results from the "worst case" pothole were used to 

predict the maximum value of A-Frame force for the howitzer. From 

this predicted maximum force the stress in the A-Frame and T-Bar 
was calculated and compared to the yield stress for the material (see 
Chapter V). 

The predicted A-Frame force is calculated by multiplying the 
combined spring stiffness of the A-Frame/T-bar times the change in 
length of A-Frame. The change in length of the A-Frame is a 

function of the difference in the angular displacements of the gun 
tube mass and the howitzer trail/saddle mass at a given time. The 
angular displacements of the gun tube and the trail were determined 

from numerically integrating the equations of motion. 

In addition to determining results in terms of A-Frame force, 

the system's mass and stiffness matrices (see Appendix B) were used 
to determine natural frequencies and mode shapes of the system. 
These natural frequencies and mode shapes were used to gain 

physical insight into the system's performance. The natural 

frequencies and mode shapes for the system are shown in Appendix 
E 
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The Matlab "FFT" function was used to perform the frequency 
analysis on the A-Frame force results from the washboard course to 
determine the frequencies at which the A-Frame force occurred in 

the model results. This frequency analysis a used 1250 Hz data 

sampling rate. This sampling rate yielded a time increment, dt, that 
equals 1/1250 or the inverse of the sampling rate. In order to have 

a frequency increment that was sufficiently small, 8192 data points 
were used for the "FFT". The record length, T, equals dt times the 
number of points or 8192/1250. The frequency increment, df, was 
the inverse of the record length or 1250/8192. 

A. WASHBOARD COURSE RESULTS 

The time histories of the computer model's predictions for A- 
Frame force from the Six Inch Washboard Course are periodic with a 
period that corresponds to the period of the washboard input. But 
unlike the sinusoidal washboard input the A-Frame force output 
predictions are not sinusoidal. These non-sinusoidal A-Frame force 

outputs most likely results from the effects of the non-linear shock 

absorber damping that was incorporated into the model (see Chapter 
II). 

The frequency analysis of the computer model’s predictions for 
A-Frame force shows that as expected the largest component of 
output force occurs at the same frequency as the washboard input 
frequency. In all cases, there is also a significant but lesser force 

component that occurs at a frequency that is twice the input 

frequency. 

The computer model was driven over the six inch washboard 

course at speeds of 3.75, 5.0 MPH, and 5.25 MPH. These speeds 
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correspond to the speeds for which actual data was taken during the 

APG test. The A-Frame force levels were relatively low at 3.75 MPH. 

The force levels in the A-Frame reached a maximum at 5.25 MPH 

because this speed for the washboard course corresponds to a 
natural frequency of the system. Driving the model at speeds higher 
than 5.25 MPH resulted in decreasing A-Frame force levels. The 
time histories and frequency analysis of the model’s predicted 
results for A-Frame force are shown in Figures 4-1 through 4-6. 

1. Results at 3.75 MPH 

The time domain model results at 3.75 MPH are shown in 

Figure 4-1. The predicted A-Frame forces are approximately 1500 

lbs in tension and 1200 lbs in compression. Figure 4-2 shows the 

frequency analysis at this speed. As expected, there is a large force 
component at the input of 3.75 MPH or a frequency of 0.92 Hz and 

there is smaller force spike at twice this frequency, 1.84 Hz. 

2. Results at 5.0 MPH 

The time domain results are shown in Figure 4-3. The A- 

Frame force levels reach 2800 lbs in tension and 4800 lbs in 

compression. Figure 4-4 shows the frequency analysis results. Again, 
a large force component at the input, 5.0 Hz, and another force spike 

at twice this frequency, 2.42 Hz. 

3. Results at 5.25 MPH 

The time domain results are shown in Figure 4-5. The A- 

Frame force levels reach 2900 lbs in tension and 4700 lbs in 

compression. Figure 4-6 shows the frequency analysis results. A 
large force component exists at the input frequency, 1.28 Hz, and 

smaller force component exists at twice this frequency, 2.56 Hz. 
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A-Frame Force 3.75 MPH 




Figure 4-1 
History) 



Model Results for A-Frame Force at 3.75 MPH (Time 



FFT at 3.75 MPH 




Figure 4-2 Model Results for A-Frame Force at 3.75 MPH (FFT) 



A-Frame Force 5.0 MPH 




Figure 4-3 Model Results for A-Frame Force at 5.0 MPH (Time 
History) 



4 6 



2000 



FFT at 5.0 M PH 




Figure 4-4 Model Results for A-Frame Force at 5.0 MPH (FFT) 



A-Frame Force 5.25 MPH 




Figure 4-5 Model Results for A-Frame Force at 5.25 MPH (Time 
History) 



FFT at 5.25 MPH 




Figure 4-6 Model Results for A-Frame Force at 5.25 MPH (FFT) 



B. NATURAL FREQUENCIES AND MODE SHAPES 

The natural, frequencies and mode shapes of the system were 

determined by using the mass and stiffness matrices that were 

developed in Chapter II. The eigenvalues and eigenvectors of the 

system are determined from the matrix that results from multiplying 
the inverse of the mass matrix, [NT *], by the stiffness matrix, [K]. 
The system's natural frequencies are the square root of the 

eigenvalues. The system's mode shapes are the eigenvectors that 
correspond to these natural frequencies. Analysis of these natural 
frequencies and mode shapes shows that the howitzer's bounce mode 
occurs at a natural frequency of 1.28 Hertz which corresponds to a 

speed of 5.25 MPH for the howitzer traveling over the washboard. 

By providing an excitation at a natural frequency of the system (1.28 
Hz) the maximum A-Frame force that could possibly be developed 
for this particular washboard input resulted. 

C. "WORST CASE" POTHOLE COURSE RESULTS 

The model results from the "worst case" pothole that was 

identified in Chapter III show that maximum A-Frame force occurs 

lbs at 15 MPH. At speeds slower or faster than this speed the model 

results showed lower levels of A-Frame force for this pothole. Figure 

4-7 shows the A-Frame force for the "worst case" pothole at 15 MPH. 

There are spikes in the A-Frame force when the wheels of the 

system hit the pothole. On Figure 4-7, at 0.1 seconds the front 

HUMMV wheels hit, at 0.6 seconds the rear HUMMV wheels hit, and 

at 1.2 seconds the howitzer wheels hit the pothole. When the 
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howitzer wheels hit at 15 MPH the maximum A-Frame force is 8,500 
lbs. 
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Figure 4-7 Model Results of A-Frame Force for "Worst Case" Pothole 
at 15 MPH 



V. STRUCTURAL ANALYSIS 



This chapter explains the structural analysis that was 

performed on the howitzer in order to determine both the force 

levels in the A-Frame and the stress on the howitzer's components. 
To determine the actual force levels in the A-Frame, the A-Frame 

was instrumented with strain gages as shown in Figure 5-1. This 

instrumented A-Frame was then calibrated to determine a 
calibration factor which would correlate strain gage readings to A- 
Frame force levels. This calibrated A-Frame was used during the 

dynamic test (see Chapter VI) at Aberdeen Proving Grounds (APG) to 

measure the actual force levels in the A-Frame. 

A. A-FRAME CALIBRATION 

The primary purpose of calibrating the A-Frame was to develop 
a reliable means to correlate dynamic strain gage measurements on 
the legs to the amount of force that causes the strain. Additionally, 

the strain measurements obtained during calibration were used to 

determine the value of the equivalent linear spring stiffness of the 
A-Frame for the computer model. During the calibration extra strain 
gages were attached to help identify levels of stress at key locations 
other than the locations needed for the force versus strain 

calibration. These extra gages provided information that was used to 
form conclusions about the stress in the A-Frame, discussed later. 




Figure 5-1 Location of A-Frame Strain Gages 




The A-Frame was loaded in both a vertical and lateral load 

configuration to determine if there were any significantly different 

effects between the two loading configurations. Figures 5-2 & 5-3 

show the test setup for both these configurations. 

1. A-Frame Instrumentation 

To conduct the calibration, the A-Frame was instrumented 

with 13 strain gages (Figure 5-1). The data taken from the pair of 

single strain gages located on the bottom of each leg was used for the 
development of the calibration curve. The pair of single strain gages 
on the top of each leg were used to verify that loads in the legs were 

mainly axial loads as opposed to bending loads. The single gage on 

the A-frame cross support brace and the single gage on the attaching 

clamp for this brace were used to identify stress levels in these 
members. The latter two strain gages were used because there were 
some known cases where the cross support brace clamps had 
cracked, and it was unknown whether these cracks were the result of 

high load levels or caused by manufacturing defects. 

2. Support Equipment 

The calibration was conducted using the lab facilities of the 
Department of Mechanical Engineering at the Naval Postgraduate 

School. To hold the A-Frame during the calibration a test support 

structure was built. The main objective of this support stand was to 

closely match the mechanical effects that occur in the A-Frame when 
it is loaded while connected to the T-Bar. The test support used was 
constructed out of heavy steel and used actual eyebolts from a Ml 19 
Howitzer to connect the bottom of the A-Frame legs to the test 
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Figure 5-2 Vertical Configuration for A-Frame Test 




Figure 5-3 Lateral Configuration for A-Frame Test 
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support (Figure 5-4). A clamping device to attach the apex of the A- 
Frame to the load cell was also designed which simulated the loading 

effects of the actual clamping device on the howitzer (Figure 5-5). 

The test support structure achieved the goal of simulating the A- 
Frame being attached to the T-Bar of the howitzer. 

3. Load Cells 

The vertical loads were applied using a Riehe Testing 
machine as shown in Figure 5-2. This loading device was used to 
apply vertical loads up to 12,000 lbs. The lateral loads were applied 
using a turn-buckle that was connected to a 5,000 pound rated 

capacity dynamometer. This dynamometer and turn-buckle system 
is shown in Figure 5-3. 

4. Calibration Procedures 

Two data runs were recorded for both the vertical and 

lateral configurations. A data run consisted of loading the A-Frame 
in increments to the maximum load for that particular data run and 
then unloading it in the same manner. All strain gages were zeroed 
before the start of each run. At each load increment, strain gage 
readings were recorded for all of the strain gages. Tables F-l 
through F-4 contain the strain gage readings and load data obtained 
from these calibration procedures. 





Figure 5-5 A-Frame Apex Clamping Device 
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Calibration Factor 



From the vertical test strain gage data in Appendix F the 
A-Frame vertical force versus strain calibration factor was 
determined. Since the A-Frame vertical load test was a static test 

the vertical force loading at the apex of the A-Frame must be 
balanced by the sum of the vertical reaction forces in the legs of the 
A-Frame (Figure 5-6). 




Figure 5-6 Balance of Forces on the A-Frame 



The reaction forces (R v i &R V 2 ) at the bottom of each leg of 

the A-Frame were determined by using the strain gage readings 
from Gage Pair 1 (gages 1 & 2) and Gage Pair 4 (gages 9 & 10) that 

are identified in Figure 5-6. The average value of the strain for 

either of the strain gage pairs was multiplied by the Modulus of 
Elasticity of the A-Frame (10 x 10 6 psi) which yields the axial stress 

in the A-Frame leg. This axial stress value was then multiplied by 
the cross sectional area of the A-Frame legs (0.6107 in 2 ) to 
determine the axial force. To calculate the vertical component of this 
axial force the above value for axial force was multiplied times the 
cosine of 28°. A sample calculation is shown below using the 10,000 
lb load increment of Vertical Test #1 (see Table F-l for strain gage 
data): 



Strain Gage Pair 1 => Avg Strain = 988 * 10‘ 6 
.-.Stress = Strain x Modulus of Elasticity 
= (988 x 10' 6 )(l0 x 10 6 psi) = 9,880 psi 
=* Force = Stress x Area = (9,880 psi)(0.6107 in 2 ) = 6,033 lbs 
.-. The vertical reaction force: R vl = (6,033 lbs)(cos(28)) = 5,327 lbs 

Strain Gage Pair 4 => Avg Strain = 1,095 x 10" 6 
.-. Stress = Strain x Modulus of Elasticity 
= (l,095 x 10' 6 )(10 x 10 6 psi) = 1,095 psi 
=»• Force = Stress x Area = (10,950 psi)(0.6107 in 2 ) = 6,687 lbs 
.-. The vertical reaction force: R v2 = (6,687 lbsXcos(28)) = 5,904 lbs 



Vertical Forces - R vl + R v2 = 5,327+5,904 = 11,231 lbs 



To determine the calibration factor for vertical loading, 

the percentage difference between the known apex load and the sum 
of the calculated vertical force components at the bottom of the A- 
Frame was determined. The average percent difference for all the 
load increments is equal to the calibration factor. This calibration 
factor can then be used to compare measured loads in the A-Frame 
legs to the applied load at the apex of the A-Frame. An average 
calibration factor of 0.88 was determined over all of the load 
increments from both vertical tests. A sample calculation of the 

calibration factor for the 10,000 lb force increment of Vertical Test 
#1 (see Table F-l) is shown below: 

Actual Vertical Force => 10,000 lbs 

Measured Vertical Force =► 11,231 lbs 

Calibration Factor => = 0.89 

11,231 

6. A-Frame Stiffness 

The measured strain data was also used to determine a n 
effective linear spring stiffness for the A-Frame. To determine this 
stiffness the applied vertical load for a specified load increment was 
divided by the corresponding amount of vertical deflection in the A- 
Frame. The average axial deflection in the A-Frame legs was 
calculated from the average strain at the bottom of each leg. The 

average axial deflection was then multiplied by the cosine of 28°. to 
calculate the amount of vertical deflection. A sample calculation for 
the 10,000 lb load increment of Vertical Test #1 (see Table F-l) is 
shown below: 
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Strain Gage Pair 1 = 998; Strain Gage Pair 4 = 1095 

Avg axial strain = = 1042 x 10" 6 ; Multiply by leg length 

=> Avg axial deflection = (28 in)(l042 x 10^) = 0.029 in 
Convert to vertical deflection => (0.029 in Xcos(28)) = 0.0257 in 

A - Frame Stiffness = = 388,400 lbs/ in 

00257 m 



7. Observations 

From analyzing the results of both the vertical and lateral 
A-Frame loadings the following conclusions were reached regarding 
A-Frame loading: 

* The strain measurements in both the cross member brace and 
the attaching clamp for this brace on the A-Frame (strain gages 3 
& 4, in Figure 5-1) were very low regardless of the force at which 
the A-Frame was loaded. This indicates that these members are 
not critical components of the A-Frame in regards to failure. 

* The A-Frame loads are carried mostly through axial loads in the 
A-Frame legs as opposed to bending loads. This means that the 
bottom connection joints where the A-Frame attaches to the T-Bar 
act as pinned joints However, because there was some bending in 
the A-frame legs at lower force levels the bottom joints do not act 
as perfect pin joints. The higher force levels break the static 
friction effects in the joints causing the joints to act more as 
perfect pin joints. 



B. STRESS ANALYSIS 

The stress analysis focused on determining the highest stress 
levels in the T-Bar and the A-Frame. To determine if mechanical 
failure will occur due to exceeding the yield stress, the stress at the 



probable points of highest stress in the A-Frame and T-Bar were 
compared to the yield stress values for the respective materials. 

1. T-Bar Stress 

The T-Bar and A-Frame assembly are shown below in 
Figure 5-7. Point A on this figure is identified as the point of highest 
stress. This is where the geometry of the T-Bar assembly changes 
and the stronger built up portion of the T-Bar ends. To simplify the 
stress analysis the T-Bar is modeled as a cantilevered beam. Figure 
5-8 below depicts the cantilevered beam modeling of the T-Bar as 
viewed from the end of the T-Bar and the rear of the howitzer. The 
point of highest stress was determined to be at the corner of the T- 
Bar at the base of the modeled cantilevered beam. One of these 
corners is identified as Point A in Figure 5-8. 




A-Frame 



High Stress 
Location 
(Point A) 



^ T-Bar 



Figure 5-7 T-Bar and A-Frame Assembly 



View from Rear of Howitzer 




Figure 5-8 T-Bar Cantilevered Beam Model 



To calculate the stress at Point A, the normal stress 

components that result when the axial force in the A-Frame leg, P, 

was resolved into components in the x, y, and z directions as shown 

on Figure 5-8. The stress at Point A is the sum of these .normal stress 
components using the principle of superposition of stress. To find 

the normal stress components first the A-Frame leg force, P, was 
projected into a force component along the axis of the T-Bar in the y 
direction, P y = sin(28)P, and a second force component in the xz 

plane, P xz = cos(28)P. The P xz component was further broken up into 
components in the x and z directions: P x = cos(27)P xz and P z = 

sin(27)P xz . To calculate the normal stress that was caused by these 

force components P X; P y , and P z well known strength of materials 
formulas for stress in a beam are applied. Therefore, applying 

these formulas the expression for the maximum stress in the T-Bar 
is: 
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a = 



(Eqn 5-1) 
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where: 

P x , P y , and P z : are defined in the above paragraph. 

L=8 inches: The distance from the base of the cantilever to the 

point of application. 

A=.9240 in 2 : The cross sectional area of the T-Bar. 

I x =.8962 in 4 : The bending moment of inertia about the x axis. 

I z =.5979 in 4 : The bending moment of inertia about the z axis. 

x =1.300 in: The distance in the x direction from the geometric 

center of the T-Bar to the outer edge of the T-Bar. 

z =1.005 in: The distance in the z direction from the geometric 

center of the T-Bar to the outer edge of the T-Bar. 

2. A-Frame Stress 

The probable point of highest stress in the A-Frame will 

occur around the 0.375 inch diameter transverse circular hole that 
positions the cross member attaching clamps. This hole is identified 
in Figure 5-7 above. At the location of this hole a stress 
concentration will exist. Figure 5-9 below contains the empirical 
formula from page 601 of [Ref. 4] that was used to determine the 

stress concentration in a pipe section with a transverse hole. 

From the above formula the stress concentration factor for 
axial loads was determined to be 3.73. Additional stress 
concentrations due to the effects of bending and torsional loads were 
neglected. The justification for neglecting the effects of bending 
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Radial hole in a hollow or tolid 18a Elastic sues*, axial tension 
circular shaft 



I 6 !* 



•d D 




For a solid shaft d = 0 



= * 7i whar 4 = A > + *•(£) + *»(£)' + *•(£)’ 




and where for d/D < 0.9 and 2 r/D < 0.45 
A, = 3.000 

A, = 2.773 4 1.529<//£) - 4.379(d/D) 2 

A* = -0.421 - 12.782d/D 4 22.781(d/D) 2 

A 4 = 16.841 -t- 16 678 d/D - 40.007(d//?) 2 



32 >/£> 



18b. Elastic suess. bending when hole o = k 

r — / n 4 _ j*\ 

is farthest from bending axis a ' 

and where for d/D < 0.9 and 2r/Z? < 0.3 
A, = 3.000 

I A : = -6.690 - 1.620 d/D 4 4.432 (d/D) 7 

A* = 44.739 + 1 0.724d/£) - 1 9.927(<//£>) 2 

1 A 4 = -53.307 - 25.998 d/D 4 43.258(d/£>) 2 



hCTC 4 = A > + A ’(£) + Aj (^) + a< (£) 



18c. EJastic stress, torsion 
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and where for d/D < 0.9 and 2 r/D < 0 4 
i A', s= 4.000 

| A 2 = —6.793 4 1.133d/Z? — 0.126{d/£>) 2 

A, = 38.382 - 7.242 d/D 4 6.495 [d/D) 7 

A 4 = -44.576 - 7.428 d/D 4 58.656<d/£>) 2 



hcc I = A, + A,(|) + *',(!)■ - *.(£)’ 



Figure 5-9 
Section 



Stress Concentration Due to a Circular Hole in a Pipe 



loads was that both the A-Frame calibration results (see Appendix F) 
and the results of the instrumented test (see Chapter VI) show that 
the joint where the A-Frame connects to the T-Bar acts as a nearly 
perfect pin joint at higher load levels. Any stress concentration 
reduction due to the spring pin that is driven into the hole were 
neglected because this stress reduction was assumed to be small. 
Therefore, maximum total stress in the A-Frame was based on only 
the axial stress component in the legs and was determined using the 
following formula: 

PK 

o = — (Eqn 5-2) 

A 

where: 

P = the axial load in the A-Frame leg. 

K = 3.73: the stress concentration factor for axial loading that is 

caused by the hole in the A-Frame leg. 

A = 0.6107 in 2 : the cross sectional area of the A-Frame leg. 
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Without this hole in the A-Frame the maximum allowable 

stress would be increased by the amount of the stress concentration 
factor, 3.73. This hole appears to only aid in assembling the A-Frame 
and serves no other function. Elimination of this hole would make 

the factor of safety more than three times greater. 

3. Worst Case Stress and Factors of Safety 

Using Eqn 5-1 and Eqn 5-2 above, the maximum stress was 
calculated for the T-Bar and the A-Frame. The maximum predicted 
A-Frame force of 8,500 lbs from driving the computer over the 
"worst case" pothole described in Chapter II was divided by two to 
calculate the force in each leg of the A-Frame. The resulting axial 
force, P, in each leg of the A-Frame equals 4,250 lbs. From Eqn 5-1 

the resulting stress in the T-Bar is 63,900 psi. From Eqn 5-2 the A- 
Frame stress equals 26,000 psi. 

The above stress levels are used to calculate the predicted 

factors of safety for the T-Bar and the A-Frame. The factor of safety 
is calculated by dividing the yield stress of the component by the 
predicted maximum load stress. The yield stress for the T-Bar and 
A-Frame are 116,000 psi and 50,000 psi, respectively. Based on the 

maximum predicted stress levels given above, the factor of safety is 
1.81 for the T-Bar and 1.92 for the A-Frame. These factors of safety 
show that both the T-Bar and the A-Frame are strong enough to 

withstand the stress from the worst case predicted loads. 
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VI. INSTRUMENTED TEST 



An instrumented dynamic test was conducted on the Ml 19 

howitzer at Aberdeen Proving Grounds (APG) in April of 1993. The 
howitzer was instrumented as described below in paragraph A. The 
dynamic testing was conducted by driving the howitzer over several 

of the APG courses at various speeds. The calibrated A-Frame 

described in Chapter V was used to measure the force levels in the 

A-Frame during the test. In most respects, the test results for A- 

Frame force are similar to the results obtained from the computer 
model. However, because certain physical effects were not included 
in the computer model, there are some differences between the 
results of the actual test and the results from the computer model. 

The results of this instrumented test confirmed that the M 1 1 9 
Howitzer could be safely transported in the firing position with the 
A-Frame attached. 

A. INSTRUMENTATION 

For this test the howitzer was instrumented with 20 strain 

gages, two tri-axial accelerometers, and two extensiometers. This 

instrumentation is explained in the following paragraphs. 

1. Accelerometers 

The howitzer was instrumented as shown in Figure 6-1. A 
tri-axial accelerometer was placed on the T-Bar to measure the 
vertical, longitudinal, and transverse accelerations of the howitzer 




Figure 6-1 Howitzer Instrumentation for APG Test 
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trail mass. This accelerometer was not placed at the center of 
gravity of the trail mass because of difficulties in attaching an 
accelerometer at this point. The location chosen on the T-Bar 
provided accelerations on the same order of magnitude and 

representative of the accelerations at the center of gravity of the 

trail. The other tri-axial accelerometer was attached on the elevating 

mass near its center of gravity. This location is shown in Figure 6-2. 

2. Single Strain Gages 

The two pair of single strain gages on the bottom of each leg 
(gage pair 1 & 4 in Figure 5-6) of the previously calibrated A-Frame 
were used to measure the strain in the A-Frame legs. The calibration 
factor determined in Chapter V was used to correlate measured 
strains in the A-Frame legs to force levels in the A-Frame. Single 

strain gages were also attached on the top and the bottom of the T- 
Bar on both sides of the howitzer. The location of these gages is 

shown in Figure 6-3. Because the highest levels of stress on the T- 
Bar are at the corners, it was not possible to place the gages at these 

points. The gages were placed as close to the high stress points as 

possible. Even though these gages could not be used to determine 

the maximum stress on the T-Bar, they could be used to roughly 
judge the level of stress. 

3. Other Gages 

Strain rosettes were placed on the bottom of the A-Frame 

legs and on the bottom howitzer trail on each side of the howitzer. 



6 9 




Figure 6-2 Location of Elevating Mass Accelerometer 
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Figure 6-3 Location of T-Bar Strain Gages 
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An extensiometer was used on each shock absorber to determine if 
the shocks were extending far enough so that the stiff rubber bump 
stops would be activated. 

B. TEST COURSES 

Several of the standard APG test courses were used for the 
instrumented test. These courses are described below. 

1. Belgian Block Course 

This course is designed using "Belgian Blocks" to produce a 
road that has high frequency random bumps. The magnitude of each 
bump is small. This course is shown in Figure 6-4. During the test 
this course was traversed at speeds ranging from 10 to 25 MPH. 




Figure 6-4 APG Belgian Block Course 



2. Six Inch Washboard Course 

This course has six inch bumps that are evenly spaced six 
feet or 72 inches apart as shown in Figure 6-5. This course produces 
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significant heave and pitch effects on wheeled vehicles even at 
relatively slow speeds. During the test this course was traversed at 
speeds ranging from three to six MPH. 




3. Cross Country Course Number 1 

This course is similar to a relatively smooth but winding 
trail road. This course was traversed at approximately 25 MPH. 

4. Cross Country Course Number 3 

This is similar to a rough trail road which has random large, 
medium, and small sized bumps. This course was traversed at 
speeds ranging from 5 to 15 MPH. 



C. TEST PROCEDURES 

Each of the above courses were traversed by the howitzer being 
towed by the HUMMV at speeds within the ranges listed. Data was 

collected on all the channels described above, using a sampling rate 

of 2,500 Hertz. Because of memory limitations in the data collection 
equipment only approximately 30 seconds worth of continuous data 
could be stored. For each course several data runs were completed 
using various speeds. For instance, on the Six Inch Washboard 

Course, six data runs were taken at speeds ranging from three to six 
MPH. 

The speeds for each data run on the Six Inch Washboard Course 
and Cross Country Course Number 3 were measured using the 
HUMMV's speedometer, because these courses were too rough to use 

more accurate devices. For the data runs, on these courses, the 

driver tried to hold the HUMMV at the specified speed while himself 
being bounced around. Because of the extremely rough road 

conditions and inaccuracies in the speedometer, the actual speeds 
ended up being different than the speeds specified in the test plan. 
Consequently, the only accurate method of determining the actual 
speed during a data run on these courses is from the video tape 
taken during the run. 

D. TEST RESULTS AND DISCUSSION 

The Belgian Block and Cross Country Course Number 1 showed 
no significant levels of strain or acceleration even at speeds up 2 5 

MPH. These tests simply proved that the howitzer could traverse 



these courses without causing damage to the howitzer. The Six Inch 

Washboard Course and Cross Country Course Number 3, however, did 

produce significant levels of strain and acceleration on the howitzer 
components. 

The test results for the Six Inch Washboard Course were 

analyzed in both the time and frequency domains. Data runs at 
speeds of 3.75, 5.0 and 5.25 MPH were analyzed for the Six Inch 

Washboard Course and compared to the predicted results from the 
computer model. The most severe data run from Cross Country 

Course Number 3 at 15 MPH was analyzed in the time domain. The 

procedures used and the analysis of these results are as follows: 

1. Time History Results 

The time history for vertical "in plane" A-Frame force in 
each case was produced by first summing the "in plane" axial forces 

from both of the A-Frame legs. The same calculation procedures as 
in Chapter V were used for determination of the axial forces in the 
legs from the strain gage data for strain gage pairs 1 and 4 (see 

Figure 5-6). The vertical force was determined by multiplying the 

axial force times the cosine of 28°. This total vertical force in the A- 
Frame legs was multiplied by the calibration factor of 0.88 that was 
determined in Chapter V. The A-Frame lateral forces were 

calculated by determining the difference in the horizontal force 

components of the axial force in each A-Frame leg. The accelerations 

were plotted directly from the measured data. 

The force and acceleration time histories for the washboard 
course (Figures 6-6 through 6-11) show that high frequency spikes 



occur on the washboard course. These spikes are caused the braking 
system being activated as the howitzer went over the bumps for all 

speeds. Also, for the 5.25 MPH data run the spade hit the tops of the 
bumps. These effects are not included in the computer model and 

result in differences between the model results and the actual 
results, particularly, for the 5.25 MPH data because the spade was 

hitting. 

Time histories for A-Frame force for speeds 3.75, 5.0 & 5.25 
MPH from the actual test (Figures 6-6, 6-8, & 6-10) are similar to the 
time histories predicted by the computer model for the same speeds 
(Figures 4-1, 4-3 & 4-5). The high frequency force spikes do not 

appear in the model's predicted results, however, the levels of the 

basic force response in the actual results closely compare to the force 
response of the model. For example, the actual results for 5.0 MPH 
(Figure 6-8) has spikes that reach a maximum force of 5,500 lbs. But 

the basic force response on average reaches approximately 2,400 lbs 

in tension and 4,300 lbs in compression. The corresponding model 
prediction (Figure 4-3) shows forces of 2,800 lbs in tension and 
4,800 lbs in compression. A similar comparison can be made at the 

3.75 and 5.25 MPH speeds. Although the actual results are not 
exactly the same as the model predictions, they are close enough to 
conclude that the model is reasonably representative of the actual 
results. 

The A-Frame vertical force results for the most severe 

portion of the Cross Country Course Number 3 are shown in Figure 6- 
12. These results show that on a rough trail road the A-Frame force 



is not as great as the force that was caused by the Six 
Washboard Course at a speed of 5.25 MPH. 



Vertical Force from APG Test (3.75 MPH) 




Figure 6-6 A-Frame Force at 3.75 MPH on the Washboard Course 
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Figure 6-7 Trail Vertical Acceleration at 3.75 MPH on the 
Washboard Course 



A-Frame Vertical Force from APG Test (5.0 MPH) 




Figure 6-8 A-Frame Force at 5.0 MPH on the Washboard Course 
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Accel (Gs) 



Vertical Trail Accel from APG Test (5.0 MPH) 
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Figure 6-9 Trail Vertical Acceleration at 5.0 MPH on the Washboard 
Course 



A-Frame Vertical Force from APG Test (5.25 MPH) 




Figure 6-10 A-Frame Force at 5.25 MPH on the Washboard Course 



Vertical Trail Accel from APG Test (5.25 MPH) 




Figure 6-11 Trail Vertical Acceleration at 5.25 MPH on the 
Washboard Course 
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Figure 6-12 A-Frame Force at 15 MPH on Cross Country Course 
Number 3 



2. Frequency Domain 

The frequency domain analysis for the actual A-Frame 
force and trail acceleration results are shown in Figures 6-13 through 
6-18. The same procedures using the "FFT" function in Matlab as 
were used in Chapter IV for the model results were used to perform 
the frequency analysis on the actual results. The time increment 
used was 1/1250 seconds and was determined by using every other 



data point 


of the 


original data that 


was 


taken at 2500 Hz. 


8192 


points were 


used 


to achieve an adequate frequency discretization. 


The record 


length 


was 8192/1250 


and 


frequency increment 


was 


1250/8192. 


The 


FFT analysis of 


both 


A-Frame force and 


trail 



vertical acceleration for each data run clearly shows two frequency 
spikes. The lowest frequency and largest spike corresponds to the 
forcing frequency or road input to the howitzer tires. The second 
largest spike corresponds to twice the forcing frequency. 

The frequency analysis for actual A-Frame force for speeds 
of 3.75, 5.0, & 5.25 MPH (Figures 6-13, 6-15, & 6-17) compares 

reasonably well to the corresponding frequency analysis from the 
model results (Figures 4-2, 4-4, & 4-6). For example, the FFT results 
from the actual data at 5.25 MPH (Figure 6-17) shows a large force 
component 1180 lbs at the road input frequency (1.28 Hz) and 
second large force component of 575 lbs at double the road input 
(2.56 Hz). Where as the FFT of the model for the corresponding 
predicted A-frame force (Figure 4-6) shows a large force component 
of 1375 lbs at the road input frequency (1.28 Hz) and a second large 
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force component of 500 lbs at twice the input frequency (2.56 Hz). 
For this case, the total amount of force from the 1.28 and 2.56 Hz 
spikes for the actual results and the model results are approximately 
the same (1765 lbs versus 1875 lbs). Although FFT results for the 
actual data and the model predictions are not exactly the same, they 
are close enough further to validate the model's predicted results. 



Vertical Force FFT from APG Test (3.75 MPH) 




Figure 6-13 FFT of A-Frame Force at 3.75 MPH on the Washboard 
Course 
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Vertical Trail Accel FFT from APG Test (3.75 MPH) 




Figure 6-14 FFT of Trail Vertical Acceleration at 3.75 MPH on the 
Washboard Course 
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Vertical Force FFT from APG Test (5.0 MPF1) 




Figure 6-15 FFT of A-Frame Force at 5.00 MPH on the Washboard 
Course 



Accel (Gs) 



Vertical Trail Accel FFT from APG Test (5.0 MPH) 




Figure 6-16 FFT of Trail Vertical Acceleration at 5.0 MPH on the 
Washboard Course 



Vertical Force FFT from APG Test(5.25 MPH) 




Figure 6-17 FFT of A-Frame Force at 5.25 MPH on the Washboard 
Course 
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Vertical Trail Accel FFT from APG Test (5.25 MPH) 




Figure 6-18 FFT of Trail Vertical Acceleration at 5.25 MPH on the 
Washboard Course 
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3. Other Results 

The strain reading from the strain rosettes that were 

attached to the bottom the howitzer's trail showed low strain 

readings for all courses. These low strain readings clearly showed 
that transporting the howitzer in the firing position with the A- 

Frame attached produces insignificant levels of stress in the howitzer 
trail. 

Lateral forces in the A-Frame were also determined from 

the actual test data. Figures 6-19 through 6-21 show the lateral 

force results for 3.75, 5.0 & 5.25 MPH, respectively. The maximum 
lateral forces were less than ten percent of the maximum vertical 
forces for all runs on the washboard course. For the rough terrain 
Cross Country Course Number Three the maximum lateral force 

(Figure 6-22) was twenty percent of the maximum vertical force. 

These relatively low values of lateral force help to show the validity 

of the modeling technique used, since the computer model only 

predicts the vertical forces in the A-Frame and assumes that the 

lateral forces would be comparatively low. 



A-Frame lateral Force from APG Test (3.75 MPH) 




Figure 6-19 Lateral Force from the APG Washboard Course at 3.75 
MPH 
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A-Frame Lateral Force from APG Test(5.0 MPH) 
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Figure 6-20 Lateral Force from the APG Washboard Course at 5.0 
MPH 



A-Frame Lateral Force from APG Test (5.25 MPH) 
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Figure 6-21 Lateral Force from the APG Washboard Course at 5.25 
MPH 
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A-Frame Lateral Force for Cross Coutry #3(15 MPH) 




Figure 6-22 Lateral Force from Cross Country Course Number Three 
at 15 MPH 
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VII. 



CONCLUSIONS /RECOMMENDATIONS 



A. CONCLUSIONS 

The computer model developed, produced output results that are 

reasonably representative of the actual results for the washboard 
type road input used at Aberdeen Proving Grounds. Therefore, the 
model is validated and can be used to predict "worst case" loading on 
the howitzer's critical parts. These "worst case" predictions for a 

trench type pothole show the Ml 19 Howitzer can be safely towed in 

the firing position with the A-Frame attached. 

As an additional result of this study, we believe that the T-Bar is 
the most critical component on the howitzer in regards to mechanical 
failure when towing in the firing position. The predicted factor of 

safety for the T-bar based on the "worst case" pothole terrain is 1.81. 

Due to the high cost of repair, making the T-Bar out of thicker steel 
would provide an increased factor of safely. 

The predicted factor of safety for the A-Frame, from the same 
"worst case" terrain, is 1.92. This indicates that the A-Frame is 
strong enough to withstand the worst case expected loads. However, 
the A-Frame is made out of aluminum and aluminum does not have 

a fatigue limit. Thus, eventually the A-Frame is expected to fatigue. 
The A-Frame's fatigue life could be greatly extended by eliminating 
the holes in its legs. 
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The strain reading taken from the trail rosettes during the 
instrumented test at APG show that the stress in the trail is 
insignificant during towing in the A-Frame position. 

B. RECOMMENDATIONS 

The US Army should remove the towing restrictions and allow 
the howitzer to be towed in the firing configuration with the A- 
Frame attached during tactical operations. 

The T-Bar should be made out of thicker steel to provide added 
assurance that it will not fail. 

The A-Frame should be re-configured for manufacture without 
the hole in its legs. 
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APPENDIX A 



SYSTEM PARAMETER VALUES 

The following table defines the HUMMV/Howitzer system 

parameters using an inches (in), pounds (lbs), and seconds (sec) 
system of units. Gravity (G) is approximated to be equal to 386.4 

in/sec 2 . The nomenclature for the parameters that are identified in 
this Appendix is the same nomenclature that is used throughout this 

paper. Refer to Figures 2-1 and 2-2 to identify the parameters 

defined below. The values listed in the table incorporate the effects 
of both sides for those components that are identical on the right and 
left side of the HUMMV or Howitzer. For example, the spring 
stiffness for each of the front HUMMV tires, Ki, is 1323.0 lbs/in. This 
value is multiplied by two to account for the combined effect of both 

front HUMMV front tires. The value listed in the following table for 
Kj is 2646 lbs/in. 



TABLE A-l COMPUTER MODEL PARAMETER VALUES 



Symbo 


Value/units 


How Determined/Notes 


Ki 


2646.0 Ibs/in 


Reference 1/HUMMV front tires stiffness 


k 2 


1908.0 Ibs/in 


Reference 1 /HUMMV front suspension stiffness 


Kj 


3544.0 Ibs/in 


Reference 1/ HUMMV rear tires stiffness 


K 4 


5040.0 Ibs/in 


Reference 1/ HUMMV rear suspension stiffness 


*6 


2000.0 Ibs/in 


Estimated /Howitzer tires stiffness 


k 7 


1056.0 Ibs/in 


Reference 2/Modeled Howitzer vertical 
suspension stiffness 


Ka 


200,000.0 Ibs/in 


Calculated/Combine A-Frame & T-Bar stiffness 


Bt 


Bic=171.6 lbs -sec/in 
Bit= 308.4 Ibs-sec/in 


Reference 1/Bic is for compression & Bit 
for tension. Values are 20% more than Ref 1 


b 2 


B2C=199.2 Ibs-sec/in 
Bzt= 480.0 Ibs-sec/in 


Reference 1 / B 2 C is for compression & Bzt 
for tension. Values are 20% more than Ref 1 


b 3 


B3c=4o. 0 Ibs-sec/m 
B3T=132.0 Ibs-sec/in 


Calculated from RIA data/ B 3 C is for 
compression & B3 t for tension. 


b 4 


250,000 lbs-in-sec/rad 


Estimated/ Torsional Damping in elevating clutch 


Mi 


1.294 Ibm 


Reference T7 includes mass or HUMMVfront 
wheels.shocks and susDension sorinas 


m 2 


1.294 Ibm 


Reference 1 / Includes mass of HUMMV rear 
wheels,shocks and suspension springs 


m 3 


21.22 Ibm 


Reference 1/ HUMMV Body 


m 4 


1.035 Ibm 


Reference 2/Includes howitzer wheels, 
shocks, and 30% of the road arm masses 


M 5 


9.834 Ibm 


Reference 2/ Howitzer trail & saddle mass 


Me 


5.952 Ibm 


Reference 2/ Howitzer elevating mass 


h 


52,680.0 lbs-in-sec 2 


Reference 1/ Moment of inertia about HUMMV CG 
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TABLE A-l COMPUTER MODEL PARAMETER VALUES (CONT) 



Symbo 


Value/units 


How Determined/Notes 


>2 


10,000 Ibs-in- sec 2 


Estimated/ Moment of inertia of Trail saddle 
mass about CG 


'3 


1 5,000 lbs-in-sec 2 


Estimated/ Moment of inertia of Elevating mass 
about CG 


Li 


83.8 in 


Reference 1 /HUMMV CG to front springs & shocks 


l 2 


46.7 in 


Reference 1 /HUMMV CG to rear springs & shocks 


L 3 


75.7 in 


Reference 1 /HUMMV CG to hitch connecting point 


l 7 


130.5 in 


Measured/ HUMMV front axle to HUMMV rear axle 


L 8 


294.5 in 


Measured/ HUMMV front axle to Howitzer axle 


Ln 


104.4 in 


Measured/ Howitzer axle to CG of Trail & Saddle mass 


Li 2 


135.0 in 


Measured/Howitzer axle to hitch connecting point 


L 1 3 


38.0 in 


Measured/Trunion to where A-Frame 
connects to the T-Bar 


L 16 


110.3 in 


Measured/Trunion to hitch connecting point 


La 


45.0 in 


Measured/ Trunion to where A-Frame 
attaches to the Elevating mass/Gun Tube 


l cg 


32.5 in 


Provide by RIA/ Tnnion to CG of Elevating mass 


00 


1 2.0 degrees 


Measured/ Static angle from horizontal to 
Gun Tube when A-Frame is A tithed 


*A 


51.0 degrees 


Measured/ Angle from Gun Tube to A-Frame 
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APPENDIX B 



MATRICES FOR EQUATIONS OF MOTION 
This Appendix contains the mass, damping, and stiffness 
matrices for the equations of motion that were developed in Chapter 
II. The mass system matrix, [M], is as follows: 



I’m,, 


0 


0 


0 


0 


0 


0 


1 


0 


m 22 


0 


0 


0 


0 


0 




0 


0 


m 33 


0 


m 35 


m 36 


m 37 




0 


0 


0 


m 44 


0 


0 


0 




0 


0 


m 53 


0 


m 5S 


m 56 


m 57 




0 


0 


m 63 


0 


m 65 


m 66 


m 67 




0 


0 


m 73 


0 


m 75 


m 76 


m 77 





where: 



j m 22 = M 2 


m 33 = M 3 + M 5 + M 6 


m 35 ~~ ^5^3 + ^6^3 


m 36“^5^11 + ^6^16 


m 37 = M 6 L cg cos( 6 0 ) 


m 44 = M 4 


m 53 = ^ 5^3 + M 6 L 3 


m 55 = ^1 + M 5 L 3 + M 6 L 3 


^56 ” 1 + ^ 6 L 3 L 16 


^57 “ ^ 6 ^CG ^3 COS(0 o ) 


m 63 = M 5 L 1 1 + M 6 L 16 


m 65 ~ M 5 L 3 Ln + M 6 L 3 Lj 6 


m 66 ~ ^2 + ^5^11 + ^6^16 


m 67 M 6 L cg L 16 COS( 0 o ) 


m 73 “ ^ 6 ^CG COS(0 o ) 


ITI 75 ” ^ 6 ^CG ^3 C °S(0 O ) 


^76 “ ^^6 ^CG^l 6 COS( 0 0 ) 


1 TI 77 — I 3 + M 6 L cg 



Note due 


to symmetry 


of [M] the 


following matrix elements are 


equal: 




m 53 m 3S 


m 63 m 36 


m 7} = m 37 




m 65 = m 56 


m 75 = m 57 


m 76 = m 67 
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The damping matrix for the system is dependent on whether 
the shock absorber forces are tensile or compressive at a given time. 
The tensile or compressive state for each of the shock absorbers is 
determined by the relative velocities of the end connection points for 



the 


shock 


absorber. 


The 


damping matrix 


rc n 


0 


Ci 3 


0 


C 1S 


0 


0 1 


0 


C 22 


C 23 


0 


C 25 


0 


0 




C 31 


C 32 


C33 


C 34 


C35 


c 36 


0 




0 


0 


C43 


C 44 


C45 


c 4 6 


0 




C 51 


C52 


C53 


C 54 


C55 


C56 


0 




0 


0 


C 63 


C 64 


c 6 s 


C 6 6 


C 67 




0 


0 


0 


0 


0 


C 76 


C 77_ 





where: 



C 11 ®1 




a 

11 

l 

W 




Ci 5 = BjLj 


C 22 = ®2 




c 23 B 2 




C 25 = ”®2^2 


O 

W 

II 

1 

Jtf 




C 32 ” “®2 




C 33 ” + ®2 + ®3 


C 34 - - B 3 




C 35 + B 2 L 2 + B 3 L 3 




C 36 = ®3^12 




C 43 ®3 




c 44 “ b 3 


C 45 = “®3^3 




C 46 “ ”®3^12 




C 51 ~ 


C 52 = “®2^2 




C 53 + ®2^2 + ^ 3^3 




C 54 = “BjLj 




Css = B.L , 2 + 


B 2 L 2 + b 3 l, 




c 56 = 




C 63 = BjL 12 




C 64 “ ”®3^i2 


C 65 = 




C 66 = ^3^12 + 


b 4 


C67 “ ~B 4 


C 76 — — B 4 




C 77 “ ®4 






Note due to 


symmetry 


of [C(x)J the 


following elements are equal: 


C 31 “ C 13 


C 32 = C 23 


C 51 “ C 15 


C 43 c 34 




C 52 = C 25 


C 53 ~ C 35 


C 54 = C 45 


C 63 ~ C 36 




C 64 = C 46 


C 65 = C 56 


C 76 = C 67 
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The system stiffness matrix, [K], is as follows: 



rk„ 


0 


^13 


0 


k 15 


0 


lo 


k 22 


k 23 


0 


k 25 


0 


! ksi 


k 32 


k 33 


k 34 


k 35 


k 36 


0 


0 


k 43 


k 44 


k 45 


k 46 


1 ksi 


^52 


k 53 


k 54 


k 55 


k 56 


|0 


0 


k 63 


k64 


k 65 


k66 


Lo 


0 


0 


0 


0 


k 76 



0 

0 

0 

0 

0 

k 



67 

k 



77 



where: 



k „=K, + K 2 


u» 

II 

1 

to 




k is K 2 L[ 


tO 

II 

* 

+ 

* 

■u 


^23 ^4 




k l 5 ^4^2 


;T 

ii 

i 

W 

to 


k 32 - -K 4 




k 33 = K 2 + K 4 +K 7 


1 

II 

rr 


^35 ^ 2 ^i 


+ k 4 l 2 + 


K 7 L 3 


k^ 6 K 7 L 12 


1 

II 

c*i 

rf 




k 44 = K 6 + K ? 
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APPENDIX C 



NUMERICAL INTEGRATION COMPUTER CODE 



This Appendix contains the computer code used to integrate 
the equations of motion using the ODE 45 integration routine in 
Matlab.. 

1. THE ODE 45 INTEGRATION FUNCTION 

function xdot=mateqn(t,x) 

%Define gravity 
G=386.4; 

%Define HUMMV Parameters: 

Ll=83.8; L2=46.7; L3=75.7; L7=130.5; L8=294.5; 

Kl=2*1323.0; K2=2*954.0; K3=2*1777.0; K4=2*2520.0; 

K5=1000000.0; 

B1C=2*1. 2*71.5; B1T=2*1.2*128.25; B2C=2*1. 2*83.0; 

B2T=2*1. 2*200.0; 

Wl=2*250.0; W2=2*250.0; W3=8200.0; 

M1=W1/G; M2=W2/G; M3=W3/G; 

11=52680; RW=5400; FW=3800; 

% Define Gun Parameters: 

L4=123.4; L5=11.6; LA=45.0; LS=42.6; LCG=32.5; 

LI 1=104.4; L12=135.0; L13=38.0; L14=17.0; L15=9.70; 

LI 6=1 10.3; 

K6=2* 1000.0; K7=2*528; KA=200000; KS=300; 

B3C=2*20.0; B3T=2*66.0; B4=.5*250000.0; 

W4=2*200; W5=1500; W6=2300; M4=W4/G; M5=W5/G; 

M6=W6/G; 

12=10000; 13=15000; 

Q0= 1 2.0*pi/ 1 80; PHIA=51.0*pi/180; PHIS=14.0*pi/180; 

AW=3800; TB=88000; % Define Mass matrice 

% Define the mass matrix 
M=[zeros(7,7)]; 
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M(1,1)=M1; M(2,2)=M2; M(3,3)=M3+M5+M6; M(3,5)=M5*L3+M6*L3; 

M(3,6)=M5*L1 1 + M6*L16; M(3,7)=M6*LCG*cos(Q0); 

M(4,4)=M4; M(5,3)=M5*L3+M6*L3; M(5,5)=I1 +M5*L3 A 2+M6 A 2; 

M(5,6)=M5*L3*L1 1+M6*L3*L16; M(5,7)=M6*L3*LCG*cos(Q0); 

M(6,3)=M5*L1 1+M6*L16; M(6,5)=M5*L3*L1 1 +M6*L3*L16; 

M(6,6)=I2+M5*L11 A 2+M6*L16 A 2+M6*L14 A 2; 
M(6,7)=M6*L16*LCG*cos(Q0); 

M(7,3)=M6*LCG*cos(Q0); M(7,5)=M6*L3*LCG*cos(Q0); 

M(7,6)=M6*L16*LCG*cos(Q0); 

M(7,7)=I3+M6*LCG A 2; 



% Define Stiffness Matrix 
K=[zeros(7,7)]; 

K(1,1)=K1+K2; K(l,3)=-K2; K(1,5)=K2*L1; 

K(2,2)=K3+K4; K(2,3)=-K4; K(2,5)=-K4*L2; 

K(3,l)=-K2; K(3,2)=-K4; K(3,3)=K2+K4+K7; K(3,4)=-K7; K(3,5)=- 

K2*L1+K4*L2+K7*L3; K(3,6)=K7*L12; 

K(4,3)=-K7; K(4,4)=K6+K7; K(4,5)=-K7*L3; K(4,6)=-K7*L12; 

K(5,1)=K2*L1; K(5,2)=-K4*L2; K(5,3)=-K2*L1+K4*L2+K7*L3; K(5,4)=- 

K7*L3; 

K(5,5)=K2*L1 A 2+K4*L2 A 2+K7*L3 A 2; K(5,6)=K7*L3*L12; 

K(6,3)=K7*L12; K(6,4)=-K7*L12; K(6,5)=K7*L3*L12; 

K(6,6)=K7*L12 a 2+KA*LA*L13*cos(Q0); 

K(6,7)=-KA*L13*LA*cos(Q0); 

K(7,6)=-KA*L13*LA*cos(Q0); 

K(7,7)=+KA*L13*LA*cos(Q0); 



%Define and velocities: 

VHF=x(10) - (Ll*x(12)); VHB=x(10) + (L2*x(12)); 

VGB=x(10) + (L3*x(12)) + (L12*x(12)); 

% Define Road Inputs: 

VMPH=5.0; VEL=VMPH* 1 7.6; PD=72/VEL; 
Ul=3*(sin(2*pi*t/PD)); 

if t < L7/VEL, U2=0;else, T2=t-(L7/VEL); 
U2=3*(sin(2*pi*T2/PD)); end 

if t < L8/VEL, U3=0;else, T3=t-(L8/VEL); 
U3=3*(sin(2*pi*T3/PD)); end 
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% Determine the shock absorber damping constants based on 
compression or tension 



% HUMMV Front shock 

if (VHF - x(8)) > 0.0, B1=B1T; else, B1=B1C; end 

% HUMMV Rear shock 

if (VHB - x(9)) > 0.0, B2=B2T; else, B2=B2C; end 
% GUN shock 

if (VGB - x(ll)) > 0.0, B3=B3T; else, B3=B3C; end 

% Define Damping Matrix 
B=[zeros(7,7)]; 

B(1,1)=B1; B(1,3)=-B1; B(1,5)=B1*L1; 

B(2,2)=B2; B(2,3)=-B2; B(2,5)=-B2*L2; 

B(3,1)=-B1; B(3,2)=-B2; B(3,3)=B1+B2+B3; B(3,4)=-B3; B(3,5)=- 

B1*L1+B2*L2+B3*L3; B(3,6)=B3*L12; 

B(4,3)=-B3; B(4,4)=B3; B(4,5)=-B3*L3; B(4,6)=-B3*L12; 

B(5,1)=B1*L1; B(5,2)=-B2*L2; B(5,3)=-B1*L1+B2*L2+B3*L3; 

B(5,4)=-B3*L3; B(5,5)=B1*L1 A 2+B2*L2 A 2+B3*L3 A 2; 

B(5,6)=B3*L3*L12; 

B(6,3)=B3*L12; B(6,4)=-B3*L12; B(6,5)=B3*L3*L12; 

B(6,6)=B3*L12 A 2+B4; 

B(6,7)=-B4; B(7,6)=-B4; B(7,7)=B4; 

%Define Forcing Vector 
F=[zeros(7,l)]; 

F(1,1)=K1*U1; F(2,1)=K3*U2; F(4,1)=K6*U3; 

% CHECK IF FRONT HUMMV WHEEL LOSES CONTACT 
if K1*(U1 - x(l)) < -FW, F(2,1)=-FW; K(1,1)=K2; end 

% CHECK IF REAR HUMMV WHEEL LOSES CONTACT 
if K3*(U2 - x(2)) < -RW, F(2,1)=-RW; K(2,2)=K4; end 

% CHECK IF HOWITZER WHEEL LOSES CONTACT 
if K6*(U3 - x(4)) < -AW, F(4,1)=-AW; K(4,4)=K7; end 

% transform into system of 1st order ode's 

id=[l zeros(l,6);0 1 zeros(l,5);0 0 1 zeros(l,4);zeros(l,3) 1 zeros(l,3)]; 
id=[id;zeros(l,4) 1 0 0;zeros(l,5) 1 0;zeros(l,6) 1]; 
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a=[zeros(7,7) id; -inv(M)*K -inv(M)*B]; 
fl=inv(M)* F; 
f=[[zeros(7,l)]; flj; 

%integrate 
xdot=a*x + f; 



2. THE DRIVER CODE FOR ODE 45 

% Set the Initial Condition 
x01=[zeros(14,l)]; 

% Integrate in one second increments 
t01=0.00; tf 1=0-999; %x01=xn(length(xn),:); 
[tl,x I]=ode45('mateqn’,t01,tfl ,x01,le-3); 
t02=1.0001; tf2=l .999; x02=xl(length(xl),;); 
[t2,x2]=ode45(’mateqn\t02,tf2,x02, le-3); 
t03=2.0001; tf3=2.999; x03=x2(length(x2),:); 

[t3,x3]=ode45('mateqn',t03,tf3,x03,le-3); 
t04=3.0001; tf4=3.999; x04=x3(length(x3),:); 

[t4,x4]=ode45('mateqn',t04,tf4,x04,le-3); 
t05=4.0001; tf5=4.999; x05=x4(length(x4),:); 

[t5,x5]=ode45('mateqn',t05,tf5,x05,le-3); 
t06=5.0001; tf6=5.999; x06=x5(length(x5),;); 
[t6,x6]=ode45(’mateqn’,t06,tf6,x06,le-3); 
%t07=6.0; tf7=6.999; x07=x6(length(x6),:); 

% [t7,x7]=ode45('mateqn’,t07,tf7,x07, le-3); 
%t08=7.0; tf8=7.999; x08=x7(length(x7),:); 

% [t8,x8]=ode45('mateqn',t08,tf8,x08,le-3); 

% Combine the solution state matrix 
xn=[x I;x2;x3;x4;x5;x6];%x7;x8l; 
tn=[tl;t2;t3;t4;t5;t6];%t7;t8]; 
save ghl xn tn 
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APPENDIX D 



POTHOLE COURSE COMPUTER CODE 

This Appendix contains the computer code that was developed 
to model the terrain inputs for a road with a trench type pothole. This 
program calculates and plots a trace of a pothole profile and the axle 
path followed by circular wheel rolling through this type pothole. This 
program works for any pothole that can be model with a flat bottom, 
and sides that can be approximated as straight lines. The portion of 
the code shown below only works for potholes that have an entrance 
height that is equal to the exit height. Additional code to determine 

the terrain inputs for a any type of pothole is on file at the Naval 
Postgraduate School. Computer code was also developed to model 
bump inputs, such a log in a road and maintained on file at the Naval 
Postgraduate School. 

COMPUTER CODE: 

function [x, y, y a, u]=pothole(xent, depth, theta 1 ,xbot,theta2, farht,x trail); 

% The call statement for this program/function is as %follows: 

% [x,y,ya,u]=pothole(xent, depth, theta l,xbot,theta2,farht,xtrail) 

% Input parameters used in the call are defined as follows: 

% depth - the vertical distance from the bottom to the level 
% entrance surface. 

% theta 1 - the angle of the decline on the entrance side, measured 

% from an extension of the level entrance surface to the decline. 

% xbot - the horizontal distance across the bottom. 

% theta2 - the angle decline on the exit side, measured like theta 1. 

% farht - the vertical distance from the bottom to the level exit 

% surface. 

% xent - the horizontal length of the level entrance region. 

% xtrail - the horizontal length of the level exit region. 

% Output parameters 

% x - a column vector containing a discrete amount of evenly spaced 'x' 
% values that were used to calculate the profile and axle path. 

% y - a column vector that contains the *y' values for the pothole 
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% profile. The bottom of the pothole is defined as "zero" y value. 

% ya - a column vector containing the axle path 'y' values. 

% u - a column vector containing the height of the bottom of the 

% wheel, 'u' is the input for the simulation model 

% Other parameters assigned in the program: 

% del tax - the V intraval spacing. 

% xent - the horizontal length of the level entrance region. 

% xtrail - the horizontal length of the level exit region. 

% rad - the radius of the wheel that is rolling through the pothole. 

% Varibles that are calculated based on the above assigned 
% parameters: 

% 

% rhtetal & rtheta2 - the theta 1 & theta2 in radians. 

% xnear & xfar - the x distance differance from the top to the 

% bottom for the near & far sides. 

% xcritl & xcrit2 - the x value where the wheel can no longer 
% follow a circular path as it goes around the near & far 

% comers(near into hole, far out of hole). 

% xcrl - the x value where the wheel hits the far side or corner 

% xcr2 - the x value from where the wheel hits the the exit region 

% from the near side or near comer. 

% xchl - the x value at the bottom comer on the entrance side. 

% xch2 - the x value at the bottom comer on the exit side. 

% xch3 - the x value at the top comer on the exit side. 

% xend - the ending x value of the pothole profile. 

% ycornl & ycom2 - the y values when the axle is following a 
% circular path around the near & far comers respectively 

%********assignment of pothole & wheel parameters****************** 
deltax=.01; rad=17.0; 



< y 0 * * * * * * * * * calculating other vsribles used ************************ 



rtheta 1 =theta 1 *pi/ 180;, 

xnear=depth/ tan(rtheta 1 ) ; , 

xchl=xnear+xent;, 

xch3=xch2+xfar;, 

xcritl =xent+rad*sin(rthetal);, 

xcrit3=xent+rad*sin(rtheta2);, 



rtheta2=theta2 *pi/ 1 80; 
xfar=farht/tan(rtheta2); 
xch2=xchl+xbot; 
xend=xch3+xtrail; 

xcrit2=xch3-rad*sin(rtheta2); 

xcrit4=xch3-rad*sin(rtheta2); 



% converting/rounding x-switch values to integer increments of 
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% deltax; keeps the output x, y,&ya vectors the 
num=xent/deltax;, num=round(num);, 
numl=xchl /deltax;, numl=round(numl);, 
num2=xch2/deltax;,num2=round(num2);, 
num3=xch3/deltax;,num3=round(num3);, 
num4=xend/deltax;,num4=round(num4);, 
xdl=xcritl-xent;,nl=xdl/deltax;,nl=round(nl);, 
n 1 *deltax; 

xd2=xch3-xcrit2;,n2=xd2/deltax;,n2=round(n2);, 

n2*deltax; 



same size 
xent=num*deltax; 
xchl=numl*deltax; 
xch2=num2*deltax; 
xch3=num3*deltax; 
xend=num4*deltax; 
xcritl=xent + 

xcrit2=xch3 - 



**************** pQ'J'J-J OLE PROFILE *************************** 



% ****** creating x vectors for pothole regions ***************** 
xl=0:deltax:xent;, sizl=size(xl); 

x2=xent+deltax:deltax:xch 1 

x3=xchl+deltax:deltax:xch2;, siz3=size(x3); 

x4=xch2 + deltax:deltax:xch3;, 

x5=xch3+deltax:deltax:xend;, siz5=size(x5); 

x=[xl x2 x3 x4 x5]; 

% ****** calculating y values for the pothole profile **************** 

yl=depth*ones(l :sizl(2)); 

y2= -tan(rthetal)*x2 +tan(rthetal)*xent + depth; 
y3=zeros( 1 :siz3(2)); 

y4= tan(rtheta2)*x4 - tan(rtheta2)*xch2 ; 
y5= farht*ones(l;siz5(2)); 
y=[yl y2 y3 y4 y5]; 



<y G * * * * * * * * * * * * ** * * j/^XLJE PATH ************************************** 



%****** ca l cu l at i n g ^jg t race j n the entrance region******** 
xal=0:deltax:xent;,sizal=size(xal); 
yal=(depth+rad)*ones(l:sizal(2)); 
if xcrit2 <= xcritl, 

%*** c hgcki n g if the axle motion is only circular***** 
if xcrit2 > xent, 

% finding the x intercept of the near & far side circular motion 
for xint=xcrit2:deltax:xch3, 

yl=depth + sqrt((rad. A 2) - ((xint-xent). A 2)); 
y2=farht + sqrt((rad. A 2) - ((xint-xch3). A 2)); 
if yl<=y2; break, end 
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else 



for xint=xent+deltax:deltax:xch3, 

yl=depth + sqrt((rad. A 2) - ((xint-xent). A 2)); 
y2=farht + sqrt((rad. A 2) - ((xint-xch3). A 2)); 
if yl<=y2; break, end 
end 
end 

% calculating the near and far side y values for the axle path 
xa2=xent+deltax:deltax:xint; 
ya2=depth + sqrt((rad. A 2) -((xa2-xent). A 2)); 
xa3=xint+deltax:deltax:xch3; 
ya3=farht + sqrt((rad. A 2) -((xa3-xch3). A 2)); 
yac=[ya2 ya3] ; 

xa4=xch3+deltax:deltax:xend;, siza4=size(xa4); 
ya4=(farht+rad)*ones(l:siza4(2)); 
ya=[yal yac ya4]; 

% checking & calculating axle path when its not just circular 

%****intializes a counter & the vector yac**** 
n=0;, nnum=(xch3-xent)/deltax;, yac=zeros(l,nnum); 

%***** C he C ks & determination of axle path**** 
for xa=xent+deltax:deltax:xch3, 
if xa <= xcritl, 

ycornl=depth + sqrt((rad. A 2) - ((xa-xent). A 2)); 

yfar=tan(rtheta2)*xa-tan(rtheta2)*xch2+rad/cos(rtheta2); 
if ycoml > rad & ycoml > yfar, yacl=ycoml; 
elseif yfar > rad, yacl=yfar; 
else yacl=rad; end 
elseif xa >= xcrit2, 

ycorn2=farht + sqrt((rad. A 2) - ((xa-xch3). A 2)); 
ynear=-tan(rtheta 1 ) *xa+ tan(rtheta 1 ) *xent. .. 
+depth+rad/cos(rtheta 1 ) ; 

if ycom2 > rad & ycom2 > ynear, yacl=ycorn2; 
elseif ynear > rad, yacl=ynear; 
else yacl=rad; end 

else 

ynear=-tan(rtheta 1 )*xa+tan(rtheta 1 ) *xent. .. 

+depth +rad/cos(rtheta 1 ) ; 

yfar=tan(rtheta2)*xa-tan(rtheta2)*xch2+rad/cos(rtheta2); 
if ynear > yfar & ynear > rad, yacl=ynear; 
elseif yfar > rad, yacl=yfar; 



else yacl=rad; end 



end 

n=n+l;, yac(l,n)=yacl; 

end 

5(c 5|C 3jc 5$< 5(C 3fC 5$C sfc 5$C S$C 5}c 5$< y VfllllCS in CXlt region 

xa4=xch3+deltax:deltax:xend;, siza4=size(xa4); 
ya4=(farht+rad)*ones( 1 :siza4(2)); 
ya=[yal yac ya4]; 

end 

% ********** PLOT OF POTHOLE PROFILE & AXLE PATH ************* 

ax=[0,xend,-50,yplot]; 

y=y-depth; ya=ya-depth;u=ya-17; 

x=x';y=y';ya=ya';u=u'; 

matrix=fx y ya u]; 

save road2dat.m matrix 

axis(ax); 

plot(x,y); hold on ; plot(x,ya) , pause, hold off 
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APPENDIX E 



NATURAL FREQUENCIES AND MODE SHAPES 

Table E-l below contains the natural frequencies and mode 
shapes of the HUMMV/Howitzer system. They were calculated using 
the procedures decsribed in Chapter IV. In Table E-l the system's 

natural frequencies are listed in the top row. The positional 

coordinate values of the mode shape which corresponds to a 

particular natural frequency are listed in the column below the 
natural frequency. The natural frequencies shown in Table E-l have 
units of hertz. The positional coordinates have units of inches for Zi, 
Z?,Z 3 , and Z 4 , and radians for the rotational cordinates 0i, 02 , and © 3 . 



TABLE E-l NATURAL FREQUENCIES AND MODE SHAPES 

Natural Frequencies 



ates 


1.286 


1.870 


3.328 


Zi 


-.0001 


.0783 


-.2392 


z 2 


.0130 


.1293 


.1425 


Z 3 


.0140 


.2057 


-.0408 


Z 4 


-.1064 


.0168 


.0202 


01 


.0002 


.0003 


.0059 


02 


-.0025 


-.0013 


-.0026 


0 3 


-.0025 


-.0013 


-.0027 



12.26 


13.61 


18.51 


23.80 


-.0212 


1.217 


.0130 


.0342 


.0076 


.0364 


- 1.222 


-.1221 


.0018 


-.0221 


.0138 


.0277 


- 1.385 


- 0.221 


-.0097 


.0341 


.0001 


.0011 


.0005 


.0025 


.0001 


-.0003 


-.0001 


-.0036 


.0003 


-.0006 


-.0011 


.0063 
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APPENDIX F 



A-FRAME CALIBRATION DATA 

The following tables contain the load versus strain data that 
was recorded during the A-Frame calibration procedures that are 
described in Chapter V. The locations of the strain gages is shown on 
Figures 5-1 and 5-6. Tables F-l and F-2 contain the data from the 
vertical load tests and Tables F-3 and F-4 contain the data from the 
lateral load tests. The method used to calculate stress values from 
the strain gage readings is described in Chapter V. 
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TABLE F-l VERTICAL TEST #1 A-FRAME CALIBRATION DATA 




TABLE F-l VERTICAL TEST #1 A-FRAME CALIBRATION DATA (CONT) 
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TABLE F-3 LATERAL TEST #1 A-FRAME CALIBRATION DATA 
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APPENDIX G 



HOWITZER ORIENTATION PICTURES AND DIAGRAMS 

This Appendix contains pictures and diagrams that identify the 
key components of the Ml 19 Howitzer that were important for this 
study. 




Figure G-l Towing in Firing Position 
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Figure G-2 Towing in Stowed Position 




Figure G-3 A-Frame/T-Bar Support 
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Figure G-4 T-Bar (View from top Down) 
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Figure G-5 Howitzer Wheels and Suspension 
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